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ABSTRACT  OF  OBJECTIVES  AND  ACCOMPLISHMENTS: 

The  following  three  research  phases  were  planned: 

(1)  Develop  the  general  framework  of  a  cloud-moisture  model  and  perform 
sensitivity  analyses  for  a  one-dimensional  cloud  model.  Analyze  the  3DNEPH  cloud 
and  earth  radiation  budget  data  from  January  and  July  1979  for  verification  of 
the  cloud  model. 

(2)  Develop  a  physically-based  three-dimensional  cloud-moisture  model.  Develop 
parameterization  schemes  for  cloud  microphysical  processes,  including  ice  phase. 

(3)  Verification  of  the  large-scale  cloud  model  with  available  cloud,  radiation 
and  synoptic  data.  Perform  numerical  experiments  to  investigate  the  effects  of 
cloud- radiation  interactions  and  feedbacks  on  large-scale  temperature  fields,  and 
the  efforts  of  infrared  radiative  cooling  and  ice  processes  on  the  prediction  of 
cloud  cover  and  cloud  liquid/ice  water  content. 

A  time-dependent,  three-dimensional,  large-scale  cloud-moisture  model  has  been 
developed  for  the  prediction  of  cloud  cover,  cloud  liquid/ice  water  content, 
precipitation,  specific  humidity  and  temperature.  The  wind  field  is  prescribed 
using  the  results  from  a  general  circulation  model.  The  model  includes 
parameterizations  of  ice  and  water  cloud  microphysical  processes,  and  solar  and 
infrared  radiative  transfer  involving  cloud  layers.  The  predicted  cloud  and 
radiation  fields  compare  reasonably  well  with  those  analyzed  from  3DNEPH  and 
earth  radiation  budget  data.  Infrared  cooling  enhances  the  formation  of  high  and 
middle  cloud  cover  and  increases  the  liquid/ice  water  content.  The  incorporation 
of  ice  processes,  on  the  other  hand,  reduces  the  total  cloud  water  due  to  the 
gravitational  settling  of  cloud  particles.  The  large-scale  cloud  model  developed 
in  this  study  should  be  of  value  to  the  development  of  general  circulation  and 
climate  models. 
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nephanalysis  ( 3DNEPH )  and  earth  radiation  budget  (ERB)/  have  been  processed  into  formats 
appropriate  for  verification.  The  96-hour  model  simulations  of  cloud  cover/  outgoing 
longwave  radiation  (OLR)  and  net  solar  flux  at  the  top  of  the  atmosphere  (TOA)/  and  cloud 
LWC  have  been  verified  against  data  analyzed  from  3DNEPH/  ERB  and  the  Nimbus-7  Scanning 
Multi-channel  Microwave  Radiometer  (SMMR)  satellite  observations/  respectively.  The 
predicted  cloud  IWC  is  compared  to  in-situ  observations  as  well  as  to  results  from  other 
studies.  The  predicted  cloud  and  radiation  results  compare  reasonably  well  with  those 
analyzed  from  satellite  data.  The  moisture/  cloud  fields  and  earth  radiation  budget  are 
physically  linked  together  by  the  cloud  microphysical  processes  simulated  in  the  large- 
scale  cloud  model/  which  further  relates  these  cloud  and  radiation  fields  to  large-scale 
temperature  and  vertical  velocity.  The  cloud  LWC  in  the  simulation  strongly  depends  on 
the  vertical  velocity  which  is  essential  to  the  condensational  process/  while  the  simulated 
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maximum  cloud  IWC  occurs  at  **  -15°C/  which  is  the  temperature  for  the  maximum  ice  crystal 
depositional  rate  observed  by  many  cloud  microphysicists.  Numerical  experiments  are 
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The  inclusion  of  radiation  produces  a  significant  change  in  temperature/  cloud  cover  and 
total  cloud  water  content/  while  the  inclusion  of  ice-phase  processes  generates  a  substan¬ 
tial  change  only  in  total  cloud  water  content.  If  the  cloud  LWC  is  sufficiently  large  to 
support  effective  Bergeron-Findeison • s  processes/  the  total  cloud  water  content  decreases/ 
indicating  that  gravitational  settling  is  an  efficient  mechanism  in  reducing  the  cloud  IWC. 
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Section  1 


INTRODUCTION 

1.1  The  Important  Role  of  Clouds  in  the 
Earth- Atmospherer  System 

Understanding  the  role  of  clouds  in  the  earth’s  radiation  and  heat  budgets 
becomes  one  of  the  foremost  tasks  in  the  atmospheric  sciences.  Clouds  may  have  an 
important  influence  on  earth  radiation  budget  in  determining  the  earth’s  climate, 
on  atmospheric  circulation  in  contributing  to  atmospheric  energetics,  and  on  the 
atmospheric  hydrological  cycle  in  generating  precipitation  which  is  of  economic  and 
societal  importance. 

Clouds  regularly  occupy  about  50%  of  the  sky  on  a  global  scale  and  arc  one 
of  the  most  important  regulators  of  the  earth’s  radiation  budget,  i.e.,  the  balance 
between  incoming  solar  radiation  and  outgoing  longwave  radiation  at  the  top  of 
the  atmosphere  (TOA),  by  which  the  earth’s  climate  is  determined.  On  the  one 
hand,  clouds  reflect  a  large  portion  of  incoming  solar  radiation  producing  a  cool¬ 
ing  effect,  referred  as  the  albedo  effect.  On  the  other,  they  trap  the  outgoing 
longwave  radiation  emitted  from  the  earth  resulting  a  warming  effect,  referred  as 
the  greenhouse  effect.  Since  the  gre  ^'use  and  albedo  effects  are  different  in 
sign  as  well  as  magnitude,  the  existence  of  clouds  may  have  a  profound  impact 
upon  the  sensitivity  of  the  climate  to  external  perturbations  such  as  a  doubling 
of  atmospheric  carbon  dioxide.  Tn  particular,  cirrus  clouds,  which  regularly  cover 
about  20%  of  the  globe,  are  optically  thin  and  nonblack.  The  outcome  of  the 
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greenhouse- versus-albedo  effects  is  intrinsically  modified  by  the  radiative  effects  of 
nonblack  cirrus  (Liou,1986). 

Exchanges  of  latent  heat  due  to  the  formation  and  dissipation  of  clouds  di¬ 
rectly  contribute  to  the  atmospheric  energetics  which  drive  global  atmospheric 
circulations.  More  importantly  but  less  directly,  the  radiative  properties  of  clouds 
produce  changes  in  radiative  heating  and  modify  the  atmospheric  temperature 
profiles,  which  in  turn  cause  changes  in  the  atmospheric  dynamic  structure.  For 
instance,  cloud-induced  radiative  heating  substantially  alters  the  equator-to-pole 
temperature  gradient  and  results  in  changes  of  zonal  mean  winds  by  the  virtue  of 
the  thermal  wind  relationship.  Ramanathan  etal.  (1983)  demonstrated  that  cirrus 
clouds  exert  a  significant  influence  on  the  meridional  gradient  in  upper-tropospheric 
radiative  heating,  which  in  turn  strengthens  the  tropospheric  jet  stream.  The 
amount  of  this  strengthening  depends  on  the  cirrus  cloud  emissivity,  which  is  a 
function  of  cloud  ice  content. 

The  presence  of  water  and  its  phase  transitions  in  clouds  uniquely  distinguish 
clouds  as  the  most  important  component  in  the  atmospheric  hydrological  cycle. 
Clouds  serve  as  a  modulator  to  transform  water  from  vapor  to  liquid  and  solid 
phases,  and  cause  precipitation  by  the  virtue  of  gravitation.  Precipitated  water 
makes  it  possible  for  the  earth  to  support  various  kinds  of  plant  and  animal  life,  in 
particular  human  beings. 

1.2  A  Review  of  Large-Scale  Cloud  Works 

Recognizing  the  importance  of  clouds,  numerous  efforts  have  been  directed  to¬ 
ward  determining  the  effects  of  clouds  on  the  atmosphere  in  two  ways.  One  is  data 
analysis  based  on  satellite  observations  such  as  earth  radiation  budget  data  and 
cloud  nephanalysis.  The  other  is  the  inclusion  of  a  more  physically  based  cloud 
parameterization  scheme  in  atmospheric  general  circulation  models  (GCMs). 
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1.2.1  Satellite  data  analysis 


Based  on  observed  satellite  radiation  data,  Cess  (1976),  Hartman  and  Short 
(1980),  Ohring  and  Clapp  (1980),  Ohring  etal.  (1981),  and  Cess  etal.  (1982) 
examined  the  sensitivity  of  radiation  budget  to  variations  in  cloudiness.  Cess 
(1976),  based  on  both  annual  and  seasonal  climatology,  suggested  that  the  albedo 
efFect,  as  a  result  of  variation  in  cloud  amounts,  is  balanced  by  the  greenhouse 
efFcct  either  globally  or  zonally.  However,  similar  studies  carried  out  by  Ohring 
and  Clapp  (1980),  Hartman  and  Short  (1980),  and  Ohring  etal.  (1981),  using 
the  scanning  radiometer  observations  of  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  satellites,  showed  that  the  albedo  efFect  is  much  greater 
than  the  greenhouse  efFect  from  changes  in  cloud  amounts;  i.e.,  the  net  radiation 
at  TOA  is  sensitive  to  changes  in  cloud  amounts.  In  these  studies,  the  estimations 
of  cloud  amount  changes  were  indirectly  derived  from  albedo  changes  rather  than 
satellite  observational  cloud  amounts. 

In  view  of  the  uncertainty  associated  with  cloud  amounts  in  these  studies, 
Koenig  etal.  (1987)  employed  the  cloud  climatology  data  derived  from  the  U.S. 
Air  Force  three-dimensional  nephanalysis  (3DNEPH)  cloud  data  base,  as  well  as 
the  broadband  radiation  budget  data  derived  from  the  Nimbus  7  Earth  Radiation 
Budget  (ERB)  measurements,  to  undertake  a  similar  study  on  regional,  zonally 
averaged,  and  hemispheric  scales.  They  compared  the  changes  in  solar,  infrared 
(IR),  and  net  fluxes  derived  from  ERB  data  with  respect  to  the  changes  in  cloud 
amounts  deduced  from  3DNEPH  cloud  data  and  found  the  albedo  effect  dominates 
the  greenhouse  efFect  on  a  regional,  zonal  and  hemispherical  basis. 

Ramanathan  etal.  (1989)  undertook  a  different  approach  to  analyze  the  Earth 
Radiation  Budget  Experiment  (ERBE)  data,  which  were  derived  from  the  obser¬ 
vations  of  three  satellites  in  different  orbits  to  include  diurnal  variations.  A  cloud 

forcing,  the  difference  between  clear-sky  and  cloudy-sky  fluxes  at  TOA,  was  defined 
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to  assess  the  radiative  effects  of  clouds  on  climate.  Ramanathan  etal.,  based  on  the 
data,  suggested  that  the  maximum  short-wave  cloud  forcing,  with  some  of  the  values 
over  100  W/m2,  occurred  in  the  mid-  to  high-latitude  oceanic  regions  coinciding 
with  regions  of  stratus  and  the  storm  tracks  of  mid-latitude  cyclones.  The  longwave 
cloud  forcing  had  peak  values  of  50  to  100  W/m2  over  tropical  regions  and  decreased 
^.vard  the  poles.  In  the  tropics,  the  longwave  and  shortwave  forcing  terms  nearly 
cancelled  each  other  out  with  the  net  forcing  within  ±  10  W/m2.  However,  it 
was  also  suggested  that  the  longwave  forcing  (heating)  happens  primarily  in  the 
troposphere,  while  the  shortwave  forcing  (cooling)  takes  place  largely  at  the  earth’s 
surface.  Thus,  strong  vertical  radiative  heating/cooling  gradients  may  be  generated 
by  the  clouds  in  the  tropics.  In  the  global  case,  the  shortwave  cloud  forcing  (-44.5 
W/m2)  due  to  the  enhancement  of  planetary  albedo,  exceeded  in  magnitude  the 
longwave  cloud  forcing  (31.3  W/m2)  resulting  from  the  greenhouse  effect  of  clouds. 
They  concluded  that  clouds  have  a  net  cooling  effect  on  the  earth. 

In  summary,  the  above  data  analysis  based  on  satellite  observations  indicates 
that  the  solar  albedo  effect  is  larger  than  the  IR  greenhouse  effect  at  TOA  due  to 
variations  in  cloud  cover. 

1.2.2  Numerical  model  works 

Satellite  and  in  situ  observations  over  the  past  two  decades  have  greatly  increased 
our  knowledge  of  cloud  radiative  effects  on  the  climate,  but  a  real  understanding  of 
how  clouds  feed  back  to  various  processes  of  the  earth-atmosphere  system  remains 
unknown  to  us  and  requires  efforts  in  improving  the  cloud  scheme  in  numerical 
models.  The  formation  of  clouds  in  GCMs  may  be  separated  into  the  diagnostic 
/  and  prognostic  methods.  The  former  is  an  empirical  approach  in  which  cloud  cover 
and  cloud  liquid  water  content  (LWC)  are  generated  empirically  from  other  model 
variables  (see,  e.g.,  Slingo  1980).  The  latter  is  a  prognostic  approach  with  an 
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additional  water  budget  equation(s)  for  prognostic  cloud  water  content  (see,  e.g., 
Sundqvist,  1978),  providing  a  more  direct  link(s)  among  the  radiative,  dynamical, 
and  hydrological  processes  within  models. 

1.2. 2.1  Diagnostic  method 

In  early  model  works,  the  formation  of  clouds  was  very  crude  and  lacked  physical 
foundation.  An  arbitrary  fractional  amount  of  cloud  was  assigned  to  a  given  layer 
where  large-scale  condensation  occurs.  All  the  excessive  water  vapor  was  condensed 
and  immediately  precipitated  onto  the  ground.  As  a  result,  there  was  no  cloud 
LWC  or  cloud  ice  remaining  in  the  clouds.  These  two  elements  are  necessary  for 
the  calculation  of  cloud  radiative  effects.  In  climatic  perturbation  studies,  Shukla 
and  Sud  (1981)  suggested  that  cloud  cover,  through  cloud  radiation  interaction, 
can  produce  changes  in  the  hydrological  cycle,  which  in  turn  results  in  substantial 
changes  in  thermal  forcing  and  dynamical  circulations. 

Wetherald  and  Manabe  (1980)  assigned  a  cloud  amount  of  80%  to  condensation 
so  that  the  model  atmosphere  would  equilibrate  at  a  realistic  temperature.  They 
demonstrated  that  the  greenhouse  and  albedo  effects  induced  by  cloud  cover  vari¬ 
ations  compensated  for  each  other.  Therefore,  they  concluded  that  cloud  feedback 
processes  are  unimportant  to  the  sensitivity  of  climate  to  a  forcing  such  as  a  change 
in  solar  constant  or  atmospheric  carbon  dioxide.  Nevertheless,  in  their  subsequent 
studies  with  realistic  geography,  Wetherald  and  Manabe  (1986,  1988)  found  that 
cloud  feedback  processes  enhance  the  sensitivity  of  the  model  climate  to  an  increase 
in  atmospheric  carbon  dioxide.  This  conclusion  is  in  agreement  with  that  drawn  by 
Hansen  etal.  (1983).  In  addition,  the  features  of  cloud  cover  change  obtained 
in  these  two  studies,  as  well  as  in  that  of  Washington  and  Meehl  (1984),  are 
qualitatively  similar. 
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More  recently,  increasing  efforts  have  been  placed  on  improving  either  the  ra¬ 
diative  properties  of  clouds  or  cloud  parameterization  schemes  in  GCMs  to  show 
the  importance  of  cloud-radiation  interaction  to  the  simulation  of  atmospheric 
circulation.  Ramanathan  etal.  (1983)  demonstrated  that  refinement  in  radiative 
properties  of  nonblack  high  cirrus  emissivity  significantly  improves  the  model  sim¬ 
ulation  of  atmospheric  temperatures  as  well  as  zonal  mean  winds.  Substantial 
increases  in  temperature  of  about  10°  K  and  2°  K  in  the  upper  troposphere  at 
the  winter  and  summer  polar  regions,  respectively,  and  a  decrease  in  temperature 
about  2°  to  3°  K  in  the  tropical  upper  troposphere  were  found  when  the  high 
cloud  emissivities  were  arbitrarily  taken  as  values  less  than  1.  These  changes,  by 
virtue  of  the  thermal  wind  relationship,  brought  the  magnitudes  of  their  simulated 
tropospheric  jets  into  much  closer  agreement  with  those  observed  in  the  atmosphere. 
Slingo  and  Slingo  (1988)  reached  the  same  conclusion  in  their  study  of  cloud 
longwave  radiative  forcing.  Following  the  cloud  scheme  proposed  by  Geleyn  (1981) 
in  a  simple  quadratic  form,  Liou  and  Zheng  (1984)  employed  a  GCM  designed  for 
medium  range  weather  prediction  and  showed  the  importance  of  cloud  and  radiative 
processes  to  the  maintenance  of  the  Hadley  circulation. 

To  improve  the  cloud  cover  scheme  in  GCMs,  Slingo(1980,  1987)  derived  an 
empirical  cloud  parameterization  scheme  from  the  studies  of  the  GARP  (Global 
Atmospheric  Research  Program)  Atlantic  Tropical  Experiment  (GATE)  data.  She 
empirically  related  the  model  cloudiness  to  the  large-scale  model  variables  such  as 
relative  humidity,  atmospheric  stability,  convective  activity  and  temperature  lapse 
rate.  In  her  works,  the  predicted  cloudiness  was  verified  against  that  derived  from 
Nimbus  7  and  showed  an  improvement  in  forecasting  tropical  and  extra-tropical 
cloudiness.  Slingo  also  suggested  that  a  cloud  parameterization  can  be  at  least  as 
important  as  that  of  radiative  transfer. 


1.2. 2. 2  Prognostic  method 

The  prognostic  method  which  requires  the  additional  prognostic  equation(s) 
for  cloud  LWC  and/or  cloud  ice  water  content  (IWC),  can  interactively  provide 
cloud  LWC  and  IWC  to  be  used  in  the  calculation  of  cloud  optical  properties. 
However,  the  formation  of  clouds  and  their  radiative  properties  depend  on  many 
microphysical  processes  that  can  not  be  determined  in  large-scale  models.  Nor  is 
it  possible  in  a  large-scale  model  to  resolve  these  microphysical  processes  for  each 
individual  cloud  particle.  It  also  is  not  easy  to  empirically  relate  these  microphysical 
properties,  i.e.,  cloud  LWC  and  IWC,  to  GCM  variables  such  as  relative  humidity, 
wind,  temperature  and  static  stability.  We  have  to  rely  on  the  bulk  quantities  which 
represent  the  various  source  and  sink  terms  in  the  water  budget  equation(s). 

Sundqvist  (1978)  made  the  first  attempt  to  develop  a  physically  based  prognostic 
large-scale  cloud  scheme  for  use  in  a  GCM.  In  his  approach,  the  cloud  LWC 
is  a  prognostic  variable  determined  by  various  source  and  sink  terms  which  are 
parameterized  in  terms  of  associated  bulk  quantities.  Also  partial  cloudiness  is 
allowed  to  form  when  a  large-scale  relative  humidity  is  less  than  100%. 

Sundqvist  (1981)  incorporated  the  proposed  cloud  scheme  into  the  European 
Center  for  Medium  Range  Weather  Forecasts  (ECMWF)  GCM  and  carried  out 
a  5-day  integration.  The  forecast  cloud  amounts  were  compared  with  satellite 
photographs,  which  showed  that  the  cloud  patterns  are  reasonably  well-simulated; 
the  predicted  cloud  LWC  were  within  an  observational  range.  Since  then,  there  is 
an  increasing  effort  to  include  the  prognostic  cloud  LWC  into  GCMs.  For  example, 
Sundqvist’s  cloud  scheme  has  been  adopted  in  the  Goddard  Institute  for  Space 
Studies  (GISS)  GCM.  (see,  Geno,  1990). 

Mitchell  etal.  (1989)  conducted  a  study  of  cloud  feedback  in  the  UK  Meteoro¬ 
logical  Office  atmospheric  general  circulation  model  which  is  gridded  by  11  vertical 
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of  50-m  ocean  mixed  layer.  A  cloud  liquid  water  budget  equation  was  included  in 
the  model.  Furthermore,  water  cloud  was  smoothly  converted  into  ice  cloud  with 
its  ice  crystal  bulk  terminal  velocity  arbitrarily  assumed  to  be  1  m/s.  According 
to  their  results,  which  were  averaged  over  five  years  after  the  simulations  reached 
equilibrium,  the  global  surface  warming  due  to  doubled  carbon  dioxide  was  reduced 
from  about  5.2°  to  2.7°  K  as  the  prognostic  water  was  included.  When  the  cloud 
radiative  properties,  which  were  calculated  based  on  the  schemes  developed  by  Liou 
and  Wittman  (1979)  for  short-wave  albedo  and  transmittances  and  Stephens  (1978) 
for  long-wave  flux  emittance,  were  allowed  to  depend  on  the  cloud  water  content, 
the  surface  warming  dropped  further,  to  only  1.9°  K.  They  concluded  that  clouds 
exert  a  substantial  negative  feedback  on  the  climate.  This  conclusion  is  clearly 
contradictory  to  that  obtained  by  Roeckner  (1988). 

Roeckner  (1988)  performed  a  30-year  simulation  using  the  Meteorological  In¬ 
stitute/University  of  Hamburg  GCM,  with  coupled  mixed-layer  ocean  and  sea-ice 
models  to  investigate  cloud  feedback  effects  with  response  to  a  2%  increase  in  the 
solar  constant  -  roughly  equivalent  to  a  C02  doubling.  He  separated  high  (non¬ 
black)  clouds  from  low  and  middle  (black)  clouds  and  calculated  cloud  radiative 
properties  in  the  GCM.  Roeckner  concluded,  based  on  his  results,  that  the  net 
cloud  feedback  is  positive  for  the  atmosphere  and  for  the  whole  planet  but  negative 
for  the  surface.  He  indicated  that  the  positive  feedback  of  the  cloud  optical  depth 
is  largely  contributed  by  the  increase  of  high  cloud  infrared  opacity,  especially  in 
the  tropics. 

Liou  etal.  (1990)  investigated  the  sensitivity  of  cloud  microphysics  and  radiative 
transfer  parameterizations  to  large-scale  cloud  processes  using  the  3-D  global  cloud 
model  developed  by  Lee  etal.  (1990).  In  the  cloud  model,  cloud  IWC  was  a  prog¬ 
nostic  variable  based  on  the  cloud  ice  water  budget  equation.  Liou  etal.  explored 
the  means  by  which  cloud  microphysics  and  radiative  transfer  associated  with  cirrus 
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clouds  were  physically  parameterized  in  terms  of  the  mean  effective  size  ( De )  and 
ice  water  path  (IWP).  The  IWP  was  given  by  the  IWC  and  model  thickness;  the  De 
could  be  related  to  the  IWC  through  their  dependence  on  the  maximum  dimension 
of  ice  crystals.  These  two  quantities,  De  and  IWP,  were  used  to  parameterize  the 
optical  depth  of  cirrus.  With  the  prognostics  of  IWC,  the  radiative  properties  of 
cirrus  clouds,  namely  solar  reflectivity  and  absorptivity  and  IR  emissivity,  were 
interactively  calculated  based  on  the  developed  parameterizations. 

They  found  a  strong  dependence  of  cirrus  radiative  properties  on  the  De.  For 
a  given  IWP,  smaller  De  reflected  more  solar  radiation,  emitted  more  IR  radiation 
and  enhanced  net  radiative  heating/cooling  at  the  cloud  top  and  bottom  than 
larger  sizes.  It  was  also  indicated  that  cirrus  clouds  are  largely  controlled  by  IR 
heating/cooling.  They  demonstrated  in  a  24-hour  simulation  that  a  variation  in 
the  De  from  50  to  75  fim  leads  to  an  increase  of  high  cloud  cover  in  some  latitudes 
by  as  much  as  4%  and,  at  the  same  time,  a  decrease  of  middle  cloud  cover  by  3-4% 
in  latitudes  between  60°  S  and  60°  N. 

1.3  The  Uncertainty  of  Cloud  Effects  in  Model  Studies 

Predictions  of  climate  change,  natural  and/or  human-induced,  are  deeply  de¬ 
graded  by  unfaithful  simulations  of  cloud  effects  in  numerical  models.  Indeed,  it 
has  been  debated  for  years  whether  clouds  cool  or  warm  the  climate.  For  example, 
in  his  climatic  perturbation  study  in  a  GCM,  Smagorinsky(1978)  suggested  that  the 
increase  in  CO2  will  enhance  evaporation  from  the  earth’s  surface  and  increase  the 
amount  of  low  clouds,  thus  exerting  a  cooling  effect  on  climate.  On  the  other  hand, 
studies  by  Road  (1978)  and  Schneider  etal.  (1978)  indicated  that  cloud  variations 
may  have  a  positive  feedback  effect  on  the  sensitivity  of  the  global  mean  climate  due 
to  decreasing  cloud  cover  in  the  models  as  the  sea  surface  temperature  increased. 

More  recently,  two  GCMs  with  explicit  cloud  water  budgets  produced  conflicting 

results  as  they  were  applied  to  climate  sensitivity  simulations.  Mitchell  etal.  (1989) 
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found  that  changes  in  the  relative  amounts  of  ice  and  liquid  clouds  were  sufficient  to 
reduce  the  equilibrium  surface  warming  by  half,  and  a  further  reduction  of  0.8°  I< 
in  surface  warming  was  contributed  by  the  negative  cloud  optical  feedbacks.  They 
concluded  the  cloud  feedback  to  climate  change  is  negative. 

On  the  other  hand,  Roeckner  (1988)  drew  the  conclusion  that  the  net  cloud 
feedback  is  positive  for  the  whole  planet.  In  his  simulation,  low  and  middle  clouds 
decrease  at  low  latitudes  but  increase  at  high  latitudes  in  response  to  a  2%  rise  in 
the  solar  constant.  He  decomposed  the  planetary  radiation  budget  into  surface  and 
atmosphere  components  due  to  changes  in  cloud  cover  and  cloud  water,  separately 
for  high  (non-black)  clouds,  and  for  low  and  middle  (black)  clouds.  Roeckner 
found  that  the  increase  of  cloud  LWC  produced  a  substantial  surface  cooling  which 
overcompensated  for  the  cloud  cover  effect  so  that  the  net  cloud  feedback  was 
negative  at  the  surface.  However,  the  atmosphere  was  heated  substantially  by  the 
increase  of  high  cloud  LWC  and  the  associated  rise  in  IR  flux  emissivity,  especially 
in  the  tropics.  He  suggested  that  changes  of  cloud  cover  and  LWC  in  low  and  middle 
clouds  are  only  of  minor  importance  for  the  atmosphere.  The  planetary  radiation 
change  is  obviously  dominated  by  cloud  LWC  and  the  associated  IR  emissivity  for 
high  cirrus  clouds. 

In  view  the  significant  differences  in  cloud  effects  simulated  by  current  GCMs, 
Cess  etal.  (1989)  made  an  intercomparison  of  cloud  feedback  in  14  GCMs  and  found 
a  wide  range  of  disagreement  that  varied  from  modest  negative  to  strong  positive 
feedback.  Most  of  this  disagreement  arises  from  inconsistency  in  representing  cloud 
effects  among  the  models.  They  defined  climate  and  clear-sky  sensitivity  parameters 
to  assess  a  climate  change  in  the  cloudy  and  clear  sky  conditions,  respectively.  In 
the  intercomparison  of  the  sensitivity  parameters  derived  from  the  GCMs,  they 
found  that  there  are  big  differences  among  the  global  climate  sensitivity  parameters, 
but  excellent  agreement  in  the  clear-sky  sensitivity  parameter.  Their  results  clearly 
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indicated  that  the  inconsistency  in  cloud  feedback  may  be  one  of  the  largest  sources 
of  the  wide  range  of  disagreement  in  the  climatic  sensitivity  studies. 

1.4  The  Importance  of  Cloud  Liquid/Ice  Water  to 

Radiative  Transfer 

While  it  is  difficult  to  accurately  quantify  the  effects  of  clouds  on  the  climate, 
atmospheric  scientists  generally  agree  that  the  parameterizations  of  cloud  water/ice 
content  and  cloud  radiative  properties  are  essential  to  the  differences  in  the  pre¬ 
dictions  of  climate  change.  Until  recently,  clouds  in  most  GCM  studies  have  been 
accounted  for  by  cloud  cover  only,  and  the  optical  properties  of  clouds  have  had  to 
be  calculated  separately  based  on  the  prescribed  values  of  cloud  liquid/ice  water 
content.  Cloud-radiation  feedbacks  in  GCM  have  only  resulted  from  changes  in 
cloud  cover.  Thus,  the  cloud  optical  feedback,  produced  by  shortwave  albedo  or  IR 
emissivity,  has  been  missing. 

To  return  this  missing  feedback  is  the  prime  objective  to  for  including  cloud 
LWC  and  IWC  as  additional  prognostic  variables  in  models.  For  instance,  it  may 
be  expected  that  an  increase  in  cloud  liquid  water  as  a  result  of  climate  warming 
could  enhance  cloud  albedo  to  produce  a  negative  cloud  feedback.  Somerville  and 
Remer  (19S4)  employed  the  temperature  dependence  of  cloud  optical  thickness  in 
a  radiative-convective  model  and  showed  that  a  significant  negative  feedback  arises 
from  the  albedo  effect  of  cloud  liquid  water.  In  their  study,  they  made  no  distinction 
of  high  thin  cirrus  which  may  result  in  a  positive  feedback  due  to  the  nonblack 
emissivity  of  high  clouds.  Roeckner  (1988),  on  the  other  hand,  separated  high 
(non-black)  clouds  from  low  and  middle  (black)  clouds  and  found  cloud  radiative 
properties  generated  a  positive  feedback  due  to  the  increased  infrared  opacity  of 
high  cirrus  clouds.  As  he  neglected  ice  phase  in  his  study,  Roeckner’s  conclusion  was 
vulnerable.  Since  cirrus  clouds  are  dominated  by  ice  crystals,  the  ice  phase  must 
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be  included  in  the  model  physics  in  order  to  accurately  calculate  cii  infrared 
opacity. 

Liou  etal.  (1990)  proposed  a  new  link  between  the  radiative  properties  of  cirrus 
clouds  and  climate  change,  showing  how  variations  in  the  De  of  cirrus  contribute  to 
changes  in  middle  and  high  clouds.  Due  to  differences  in  cirrus  radiative  properties 
for  different  De,  high  cirrus  clouds  that  are  composed  of  smaller  De  would  produce 
more  cooling  and  heating  above  and  below  the  cloud  region,  respectively,  than 
larger  De.  They  showed  that  more  high  cloud  and  less  middle  cloud  are  simulated 
as  the  De  increased  from  50  to  75  fim.  This  may  lead  to  the  strength  of  the 
greenhouse  effect,  since  high  and  middle  clouds  are  dominated  by  the  greenhouse 
and  albedo  effects,  respectively.  However,  it  depends  on  how  the  De  responds  to 
climate  warming  to  determine  whether  the  contribution  of  De  to  cloud  feedback  on 
climate  change  is  positive  or  negative. 

Based  on  Iieymsfield  and  Platt  (1984),  ice  crystal  sizes  and  the  consequence  of 
Dc  are  strongly  dependent  on  the  temperature,  in  that  the  higher  the  temperature, 
the  larger  the  ice  crystals  and  De.  Following  this  relation,  one  may  expect  the  size 
of  De  in  cirrus  increases  in  response  to  climate  warming  and  results  in  the  decrease 
and  increase  of  high  and  middle  clouds,  respectively.  This  may  lead  to  a  reduction 
in  the  greenhouse  effect  and,  at  the  same  time,  an  enhancement  of  the  albedo  effect. 
As  a  result,  climate  warming  may  be  neutralized  due  to  the  effects  of  De  on  the 
radiative  properties  of  cirrus  clouds. 

According  to  the  preceding  literature  reviews,  it  is  clear  that  the  effects  of  clouds 
on  the  earth’s  climate  depend  not  only  on  cloud  cover  but  also  on  the  cloud  optical 
properties  determined  by  cloud  LWC  as  well  as  cloud  IWG.  In  particular,  it  is 
well-known  that  the  radiative  effects  of  cirrus  clouds,  which  consist  of  nonspherical 
ice  crystals,  are  a  potentially  important  component  of  weather  and  climate  systems 
(  Liou, 1986).  In  order  to  include  cloud  effects  faithfully  in  the  GCM  studies,  not 


12 


only  the  cloud  cover  but  cloud  LWC  and  IWC  need  to  be  explicitly  included  in  the 
model  physics. 

1.5  The  Present  Large-Scale  Cloud  Study 

The  objective  of  this  study  is  to  develop  and  to  verify  a  prognostic  large-scale 
stratiform  cloud  model.  In  addition  to  cloud  cover,  the  model  predicts  cloud  LWC 
and  IWC  to  provide  two  more  physical  elements  for  the  calculation  of  cloud  radiative 
properties,  such  as  shortwave  albedo  or  IR  flux  emissivity.  The  model,  consisting 
of  liquid  and  ice  water  substances  budget  equations,  is  based  on  the  laws  of  mass 
conservation  and  thermodynamics.  The  ice  phase  (cloud  IWC),  which  is  physically 
included  for  the  first  time  in  a  large-scale  cloud  model,  is  used  to  account  for  cirrus 
clouds.  Various  source  and  sink  terms  are  parameterized  in  terms  of  so-called  bulk 
quantities.  The  ice  crystals  bulk  terminal  velocity  is  directly  related  to  cloud  IWC 
and  is  interactive  in  the  cloud  model.  The  bulk  parameterizations  generally  follow 
Kessler  (1969),  and  Sundqvist  (1978,  1988)  for  the  liquid  phase  and  Koenig  and 
Murray  (1976),  eymsfield  and  Platt(1984),  and  Starr  and  Cox  (1985)  for  the  icc 
phase. 

An  assessment  of  the  validity  of  any  new  parameterization  scheme  is  at  least 
as  important  as  that  of  the  development  of  the  scheme.  In  addition  to  conven¬ 
tional  meteorological  elements  such  as  temperature  and  precipitation,  the  model 
predicted  cloud  cover,  cloud  LWC  and  IWC,  as  well  as  the  radiation  budget,  need 
to  be  carefully  verified  against  observations.  Cloud  cover  observational  data  can 
either  be  compiled  from  conventional  surface  observations  or  deduced  from  satellite 
remote-sensing  data.  Ground  based  cloud  cover  observations  have  been  recorded 
for  centuries.  The  cloud  cover  climatology  compiled  by  London  (1957)  based  on 
surface  observations  is  a  widely  used  cloud  climatology. 

Aware  of  the  importance  of  satellite  data  for  climatic  studies,  a  number  of  satel¬ 
lite  observational  programs  had  been  undertaken,  e.g.,  the  International  Satellite 
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Cloud  Climatology  Project  (ISCCP,  1985),  the  First  ISCCP  Regional  Experiment 
(FIRE,  1985),  and  the  Earth  Radiation  Budget  Experiment  (ERBE,  1984).  In 
addition,  earth  radiation  budget  (ERB)  data  was  recorded  by  the  Nimbus-7  satellite 
(Kyle  etal.,  1985).  The  U.S.  Air  Force  three-dimensional  nephanalysis  (3DNEPII) 
program  is  the  only  operational,  high  resolution  cloud  analysis  based  on  both 
satellite  and  conventional  cloud  observations.  The  ERB  and  3DNEPH  data  are 
used  as  verification  data  in  this  study  and  will  be  further  described  in  Section  4. 

Compared  to  cloud  cover,  data  on  cloud  LWC  are  very  limited.  The  widely  used 
cloud  LWC  data  are  those  derived  from  the  Nimbus  7  Scanning  Multi-Channel 
Microwave  Radiometer  (SMMR).  The  SMMR  data  consist  of  vertically  integrated 
liquid  water  content  in  clouds  and  rain  over  the  oceans  equatorward  of  60°  latitude. 
The  vertically  integrated  liquid  water  content  is  referred  as  the  liquid  water  path 
(LWP).  So  far,  there  is  no  observational  cloud  IWC  data  that  can  be  utilized 
as  verification  data  on  the  global  basis.  Aircraft  in  situ  measurements  of  cirrus 
cloud  properties  such  as  cloud  IWC  were  reported  by  Ileymsfield  (1975,  1977). 
These  data  may  be  used  to  qualitatively  assess  the  model-predicted  cloud  IWC.  The 
examination  of  verification  data  based  on  satellite  observation  will  be  described  in 
some  detail  in  Section  4. 
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Section  2 


MODEL  DESCRIPTION 
2.1  Governing  Equations 

We  have  developed  a  time-dependent,  three-dimensional  cloud-moisture  model 
for  the  prediction  of  water  vapor,  cloud  LWC,  cloud  IWC,  precipitation  and  tem¬ 
perature.  The  first  four  variables  are  governed  by  the  law  of  mass  conservation.  The 
latent  heat  exchanges  during  the  phase  transitions  of  water  substances  are  ruled 
by  the  first  law  of  thermodynamics  through  the  temperature  field.  The  governing 
equations  may  be  written  as  follow: 

+  ^Tx(pfa)  +  +  T^f' »> +  ^ 

=  KvQ*  +  Qr)  -  4(1  -  v)&  +  -  4(1  -  v)S,  +  Fj,  +  Ff,  (2.1) 

Op  Op  Op 

+  ^lTx{m) + 1  +  Tz 

=  -^PnQ.  +  P(1  -  n)Er  +  (1  -  v)S,  +  FS  +  F}-  ,  (2.2) 

$  I  ^  ^  ^ 

Jtm  +  !cosi[T\(m  +  W(”5c°sffl  +  al  l”(*  ■  ”,)1 

=  ^mQ'-mP-iSBF+F1j  +  F$  ,  (2.3) 
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d  1  ft  ft  ft 

Ftm  +  +  T^mms4,)  1  +  Fm<b  ~  ^t)! 

=  ^wOi  +  >)SBF-(l-,)$,  +  Fg  +  fJ‘  ,  (2.4) 

-  Jj(-P)  =  PTlP  -  (1  -  ri)pE,  ,  (2.5) 

Equation  (2.5)  is  derived  from  the  rain  water  budget  equation  based  on  the  con¬ 
sideration  of  scaling  parameters.  The  derivation  is  shown  in  Appendix  A.  In  these 
equations,  an  overbar  denotes  the  average  of  a  variable  over  a  grid  space.  The 
variables  u,  t),  and  w  are  wind  fields  and  are  obtained  from  a  GCM;  A  and  $  are 
horizontal  coordinate  in  the  longitudinal  and  latitudinal  direction,  respectively;  p 
is  the  air  density,  T  the  temperature,  rj  the  cloud  fraction,  Cv  the  specific  heat  at 
constant  pressure,  I'd  the  adiabatic  lapse  rate,  P  the  precipitation  flux,  q„,  qi,  (ji 
the  mixing  ratios  for  water  vapor,  cloud  water,  and  cloud  ice,  respectively,  and  Lc, 
L,t  and  Lj  the  latent  heat  of  condensation,  sublimation,  and  melting,  respectively. 
Also  in  these  equations,  Qc ,  Q{  and  Qr  represent  the  condensational  heating  rales 
for  cloud  water,  cloud  ice,  and  the  radiative  heating  rate,  respectively,  P  is  the 
rate  of  autoconversion  from  cloud  droplets  to  raindrops,  Er  is  the  evaporation  of 
rain  water,  Sbf  and  S,  are  the  rates  of  change  of  the  mass  of  ice  particles  due  to 
the  Bergeron- Findeison’s  processes  and  sublimation,  respectively,  and  u>;  and  Vy 
are  the  bulk  terminal  velocities  for  cloud  water  and  ice  particles,  respectively.  In 
addition  the  vertical  eddy  flux  and  horizontal  diffusion  terms  have  been  denoted  by 
I?y  and  Fj),  respectively,  where  x  can  be  T,  qv ,  qi  or 

Because  the  radiative  heating  term  (Qr)  appears  only  in  the  thermodynamic 
equation  (2.1),  radiative  heating  affects  cloud  fields  through  the  temperature  field. 
It  is  noted  that  the  precipitation  in  Equation  (  2.5  )  is  neither  advected  nor 
accumulated,  since  we  have  assumed  all  of  the  rain  water  which  is  converted  from 
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the  cloud  LWC  is  precipitated  out  of  the  grid  box  immediately.  The  variables  to 
be  solved  from  the  Equations  (2.1)-(2.5)  are  T ,  qv ,  qi ,  <?,•,  P  and  q.  Thus,  one  more 
equation  for  the  could  cover,  77,  is  needed  to  close  the  model  equations. 

The  cloud  cover,  77,  is  derived  based  on  an  average  of  the  specific  humidities  over 
clear  and  cloudy  regions  in  the  form 

9  =  99c  +  (1  -  9)90  » 

where  qc  and  q0  are,  respectively,  the  specific  humidity  in  the  cloudy  and  clear 
regions.  After  rearranging  terms,  the  cloud  cover  is  given  by 

v  _  ?»/?»  ~  9o/ft 

9c/ 9s  —  90/9* 

Here,  we  have  introduced  the  grid-averaged  saturation  specific  humidity,  qa{T), 
which  is  a  function  of  the  grid-averaged  temperature,  T.  For  a  stratiform  cloud, 
the  temperature  within  the  could  is  approximately  equal  to  the  grid-averaged 
temperature,  i.e.,  qc  =  qa(T).  Thus, 

_  ~  (k  ~  hQ) 

,  _  (1-/10)  ’ 

where  h  =  qv/qs  and  h0  =  q0/qs  are  the  relative  humidity  in  the  grid  average  and  in 
the  clear  region,  respectively.  The  h0,  threshold  relative  humidity,  is  a  parameter 
to  be  prescribed  to  close  the  model  equations. 

The  cloud  cover  in  Equation  (2.6)  shows  a  linear  relation  between  the  relative 
humidity  and  cloud  cover.  A  similar  relation  was  found  in  the  empirical  expression 
proposed  by  Smagorinsky  (1960)  for  use  in  GCM,  as  well  as  the  results  derived  by 
Krueger  etal.  (1989)  in  meso-scale  modeling  studies  (Figure  2.1).  The  values  of  the 
threshold  relative  humidity  in  Smagorinsky  (1960)  are  adopted  in  this  study. 

With  the  introduction  of  the  threshold  relative  humidity  (h0  <  1),  partial 
cloudiness  in  the  grid  box  is  allowed  when  the  large-scale  relative  humidity  is  less 
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relationship  between  cloud  cover  and  relative  humidity. 


than  100%  .  The  cloudy  area  and  clear  region  in  a  grid  box  are  denoted  by  77 
and  (1  —  77),  respectively.  Condensation  occurs  in  the  cloudy  area;  evaporation 
takes  place  in  the  clear  region.  The  raindrops  formed  in  a  given  grid  box  are  not 
subjected  to  evaporation  within  that  box.  Instead,  the  drops  that  fall  through  this 
layer  from  higher  clouds  may  evaporate  in  the  clear  region  of  the  box.  Since  the 
averaged  terminal  velocity  of  cloud  droplets  is  on  the  order  of  1  cm/s,  which  is  close 
to  the  large  scale  vertical  velocity,  it  would  be  expected  that  the  terminal  velocity 
of  cloud  droplets  and  vertical  air  motion  roughly  cancel  each  other  in  stratiform 
clouds.  Therefore,  vertical  advection  of  cloud  LWC  in  Equation  (2.3)  is  omitted  in 
the  model  calculations. 

The  vertical  eddy  terms  are  denoted  as  Fy,  where  x  can  be  T,  ft, ,  <7/  or  ft.  The 
Fy  may  be  parameterized  in  terms  of  the  Reynold  stress,  r,  as 


where  the  Reynold  stress  is  further  expressed  by 

Here  I(y  denotes  the  vertical  eddy  diffusion  coefficients  for  the  element  of  x- 
Following  Washington  and  Williamson  (1977),  Ky  =  Ky  is  assumed  and  given 
by 

'  10s  +  106{1  -  eap[1.2  x  105(f  -  7c)]}  if  g  <  7c 
Kv  =  (2.7) 

-  +  in3  if  Si  >  ~ 

where  the  countergradient  lapse  rate,  7c,  is  taken  to  be  10-5  °K/cm,  and  0  is  a 
potential  temperature  and  given  by 

0  =  f{Po/P)%  , 
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where  P  is  air  pressure  and  Pq  is  a  reference  pressure  at  1000  mb.  The  Richardson 
number,  i?,-,  in  Equation  (2.7)  is  defined  by 


ft  =  |[(f)-7j/[(|)2  +  (|)J+10-]. 

The  horizontal  diffusion  term,  is  defined  as  follows  : 

F”  =  +  lfy(pK“cos<i>lfy  ’ 

where  Kf{  is  the  horizontal  diffusivity  of  the  element  of  x •  The  horizontal  diffusivity 
is  taken  to  be  equal  to  the  horizontal  momentum  diffusion  coefficient,  Kjj. 

According  to  Washington  and  Williamson  (1977),  the  horizontal  diffusion  coef¬ 
ficient  is  given  by 

I<h  =  2 kH'D  , 

where  k0  is  a  proportionality  constant  determined  empirically,  /  is  equal  to  the 
horizontal  grid  size,  and  D  is  the  deformation  of  flow  which  is  given  by 

d  =  (  d2t  +  dD'I*  , 


where 


Dy  — 


Ds  = 


1  .du  d 

acos^>ld  A  d(f> 

1 

acosif) 


\-+- 
[dx T  ds 


(vcos<f>)] 

(ucoscj))] 
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2.2  Cloud  Microphysical  Processes  Simulated  in  the 

Model 

The  cloud  microphysical  processes  simulated  in  the  model  are  illustrated  in 
Figure  2.2.  Both  the  liquid  and  ice  phase  processes  are  included  with  the  liquid 
phase  on  the  left  hand  side  and  ice  phase  on  the  right  hand  side  of  the  diagram- 

The  liquid  phase  processes  consist  of  condensation  of  water  vapor  into  cloud  liquid 

20 


water  at  temperatures  warmer  than  —  40°  C,  conversion  of  cloud  liquid  water 
into  rain  water  by  means  of  autoconversion,  and  precipitation  of  rain  water  onto 
the  ground.  The  evaporation  of  cloud  and  rain  water  increases  the  water  vapor. 
The  parameterization  of  autoconversion  rates  follows  the  method  that  was  first 
proposed  by  Sundqvist  (1978)  to  combine  the  coalescence  and  accretion  processes 
into  a  simple  exponential  form. 

The  ice  phase  processes  include  homogeneous  nucleation  of  water  condensate 
to  cloud  ice  at  temperatures  colder  than  —  40°  C,  and  conversion  of  cloud  ice 
to  precipitating  ice  by  the  virtue  of  gravitational  settling,  which  is  a  method  of 
calculating  the  downward  ice  flux  by  integrating  tin.  mass  and  falling  velocity  of 
an  individual  ice  crystal  over  the  whole  spectrum  of  ice  crystals.  The  sublimation 
of  cloud  ice  increases  the  water  vapor.  At  temperatures  above  the  freezing  point, 
cloud  ice  melts  into  cloud  water  by  the  melting  process. 

An  important  link  between  the  liquid  and  ice  phases  is  so-called  Bergeron- 
Findeison’s  processes.  They  are  based  on  the  fact  that  the  saturation  vapor  pressure 
over  ice  is  less  than  that  over  water.  As  a  result,  the  ice  crystals  in  a  super¬ 
cooled  cloud  grow  by  vapor  diffusion  at  the  expense  of  the  supercooled  liquid  cloud 
droplets.  Once  generated  by  ice  nucleation,  ice  crystals  grow  by  deposition  due 
to  Bergeron-Findeison’s  processes  and  are  depleted  by  sublimation  as  well  as  the 
gravitational  settling.  Above  the  freezing  point,  clouds  consist  of  cloud  water  and 
contain  no  cloud  ice.  In  the  region  0°  to  —  40°  C,  cloud  water  and  cloud  ice 
may  coexist.  At  temperatures  below  —  40°  C,  cloud  ice  dominates  because  of 
homogeneous  nucleation. 

Formation  of  rimed  crystals,  such  as  graupel  and  hail,  is  not  considered  in  this 
model  since  it  occurs  in  a  small  spatial  scale.  Furthermore,  graupel  and  hailstone  arc 
observed  more  frequently  in  convective  clouds  than  in  large-scale  stratiform  clouds. 
The  parameterizations  of  source  and  sink  terms  in  the  microphysical  processes  are 
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Figure  2.2.  Cloud  microphysical  processes  simulated  in  the  model. 
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presented  at  greater  length  in  Section  3. 

2.3  Numerical  Implementation 

In  the  vertical,  a  16-layer  stretch  z-coordinate  is  used  in  the  cloud  model.  It  has 
fine  resolution  in  the  lower  levels,  where  clouds  start  to  form,  and  coarse  resolution 
in  the  upper  levels  of  the  model.  The  vertical  grid  structure  and  the  distribution 
of  variables  are  shown  in  Figure  2.3.  In  the  figure,  a  dashed  line  denotes  a  layer; 
the  interfaces  between  layers  are  referred  to  as  levels.  The  values  of  tD,  Ky,  and  P 
are  placed  at  levels  while  the  rest  of  variables  are  located  at  layers. 

In  the  horizontal,  the  grid  point  space  is  characterized  by  48  equally-spaced  lon¬ 
gitudinal  points  and  38  Gaussian  latitudes,  which  is  the  same  as  the  horizontal  grid 
structure  used  in  the  Air  Force  Geophysics  Laboratory  (AFGL)  GCM,  and  which 
provides  the  cloud  model  with  initial  data.  As  a  result,  the  data  interpolations  in 
the  horizontal  between  the  two  model  are  unnecessary.  The  horizontal  grid  space 
in  the  model  is  modified  near  the  poles  to  maintain  a  nearly  constant  geographical 
distance  between  the  grid  points.  This  is  achieved  by  decreasing  the  number  of 
points  around  the  latitude  circles  polarward  of  60°  so  that  the  distance  between 
points  is  close  to  that  at  60°  latitude. 

The  upstream  scheme,  a  first  order  approximation  in  time  and  space,  is  utilized 
for  numerically  differencing  the  horizontal  and  vertical  advection  terms.  The  scheme 
is  conditionally  stable  and  its  stability  criterion  follows  the  Courant-Friedrichs-Levy 
(CFL)  condition,  i.e.,  CmaxAt  /  AS,  where  Cmax  is  the  maximum  velocity  in  the 
computational  domain,  A t  is  the  time  step  and  AS  is  the  grid  size.  This  criterion  is 
easily  satisfied  in  large  scale  models.  For  example,  in  the  vertical,  wmax  ~  10°cm/s, 
At  ~  103  second,  and  AS  ~  lO^cm,  the  criterion  is  much  less  than  1.  The  time 
step  used  in  this  study  is  30  minutes.  The  radiation  calculations  are  updated  every 
3  hours,  i.e.,  the  radiative  heatings  are  calculated  every  3  hours  at  every  horizontal 

giid  point  and  these  values  are  used  at  each  time  step  in  the  following  3-hour  period. 
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Figure  2.3.  Cloud  model  vertical  grid  structure  in  Z-coordinate. 
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Advantages  of  the  upstream  scheme  are  to  preserve  the  positive  definiteness 
of  the  initial  values  and  produce  no  dispersive  ripples.  These  numerical  ripples 
are  inevitable  in  a  higher-order  accuracy  scheme  and  may  eventually  lead  to  the 
instability  of  the  whole  scheme.  In  addition,  the  upstream  scheme  is  less  time 
consuming  than  high-order  schemes.  However,  the  upstream  scheme  is  known  to 
suffer  from  strong  numerical  diffusion.  To  overcome  this  problem,  Smolarkiewicz 
(1983,  1984)  introduced  the  ‘antidiffusion  velocity’  into  the  basis  of  the  upstream 
scheme  and  developed  a  positive  definite  advection  scheme  with  small  numerical 
diffusion.  The  developed  scheme,  sharing  the  same  finite  difference  form  as  the 
upstream  scheme,  is  second-  or  higher-order-accurate  in  space  depending  upon  the 
implementation  of  the  antidiffusion  velocity  as  well  as  the  number  of  iterations  in 
the  corrective  steps.  We  plan  to  update  the  upstream  scheme  with  the  antidiffusion 
velocity  included  in  the  model  in  the  future. 

2.4  Radiative  Transfer  Parameterization 

The  radiative  transfer  parameterization  program  used  in  this  study  follows  that 
developed  at  the  University  of  Utah.  It  is  based  on  a  broadband  method  and 
involves  the  transfer  of  thermal  IR  and  solar  radiation  in  clear  and  cloudy  regions. 
In  a  clear  atmosphere,  the  entire  IR  spectrum  is  divided  into  five  bands:  three 
for  II2O,  one  for  CO2,  and  one  for  O3  absorption.  The  parameterizations  of  these 
broadband  IR  emissivities  were  developed  by  Liou  and  Ou(1981),  and  Ou  and 
Liou(1983).  The  solar  spectrum  consists  of  25  bands:  six  for  II2O,  one  for  CO2, 
which  overlaps  the  li20  2.7/rm  band,  and  18  for  O3.  Parameterizations  of  the 
broadband  solar  absorptivities  for  these  absorbers  were  documented  in  Liou  etal. 
(1984). 

In  a  cloudy  atmosphere,  low  and  middle  clouds  are  treated  as  blackbodies  in  the 
IR  radiative  transfer  calculation.  The  broadband  IR  emissivity,  reflectivity,  and 
transmissivity  for  high  clouds  as  well  as  the  broadband  solar  absorption,  reflection, 
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and  transmission  values  for  various  cloud  types  are  computed  based  on  the  cloud 
LWC  and  IWC  interactively  generated  by  the  cloud  model.  The  cloud  radiative 
properties  are  calculated  based  on  the  parameterizations  developed  by  Liou  and 
YVittman  (1979).  Accuracy  of  all  the  aforementioned  parameterizations  was  verified 
via  more  comprehensive  and  exact  radiative  transfer  calculations  described  in  the 
report  by  Ou  and  Liou  (1988). 

A  total  cloud  cover  is  used  to  obtain  linearly  weighted  radiative  fluxes  between 
the  clear  and  cloudy  regions.  By  interactively  providing  the  cloud  cover,  cloud 
LWC  and  IWC  to  the  radiative  transfer  calculation,  clouds  can  alter  the  radiative 
heating  profiles,  which  in  turn  feed  back  to  the  clouds  through  the  thermodynamic 
equation.  Thus,  cloud- radiation  interaction  through  the  medium  of  the  cloud  cover, 
cloud  LWC  and  IWC  is  accomplished. 

Thf*  -adiative  heating  rate  at  level  z  is  related  to  the  divergence  of  the  net  fluxes 
and  is  given  by 

dT  1  x 


where  A z  and  A F  are  the  model  layer  thickness  and  the  net  radiative  flux  difference 
between  the  layer  top  and  bottom,  respectively.  This  equation  applies  to  both 
111  and  solar  heating  rate  calculations  with  downward  solar  flux  and  upward  IR 
flux  defined  as  positive.  For  the  solar  radiation  in  a  cloud  layer,  the  net  flux 
density  decreases  from  the  the  cloud  top  to  the  bottom  because  of  cloud  absorption. 
The  absorbed  radiate  energy  heats  the  cloud  layer.  Thus,  solar  radiative  heating- 
increases  the  atmospheric  temperature.  For  a  black  cloud  (a  low  or  middle  cloud), 
the  IR  radiative  flux  is  proportional  to  the  fourth  power  of  the  absolute  temperature 
of  the  cloud.  Since  the  atmospheric  temperature  usually  decreases  with  height,  the 
downward  solar  flux  emitted  from  the  atmosphere  above  a  black  cloud  top  is  much 
smaller  than  that  emitted  upward  by  the  black  cloud.  As  a  result,  the  net  flux 
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near  the  cloud  top  is  greater  than  zero  and  leads  to  strong  IR  cooling  in  the  region. 
Similarly,  there  is  IR  heating  near  the  cloud  bottom.  For  a  nonblack  high  cloud, 
the  intracloud  flux  profile  is  determined  by  exponential  interpolation  between  the 
cloud  top  and  base  fluxes. 

Specification  of  the  overlap  that  occurs  between  clouds  layers  is  required  for  the 
calculation  of  radiative  transfer  in  a  radiation  model,  since  only  a  cloud  fraction 
for  each  model  layer  is  determined  by  the  cloud  model.  In  this  study,  the  model¬ 
generated  clouds  are  strapped  into  three  cloud  decks.  The  parameterization  of 
radiative  flux  exchanges  in  the  three-deck  cloudy  atmosphere  has  been  documented 
in  Ou  and  Liou  (1988). 

For  the  AFGL  GCM,  low,  middle,  and  high  clouds  fill,  respectively,  layers 
(4,  5,  6),  layers  (7,  8),  and  layers  (9,  10),  as  indicated  in  Figure  2.4.  Each  deck 
of  clouds  is  the  respective  averaged  cloud  cover  of  the  model  layers.  For  example, 

m  =  ^(i?4  +  Vs  +  Vg)  . 

1, 

Vm  =  ^{V7  +  Vs)  » 

Vh  =  ^{V9  +  V  10)  • 

Where  ift,  i)m  and  %  denote  deck  cloud  cover  for  low,  middle,  and  high  clouds, 
respectively. 

To  determine  the  total  cloud  amount,  we  Ihave  assumed  that  clouds  overlap 
each  other  in  a  statistically  random  manner,  which  is  consistent  with  the  model 
statistically  diagnostic  cloud  cover  scheme  (  Equation  2.6).  Thus  the  total  cloud 
cover  may  be  given  by 

T]t  =  1  -  (1  -  ry)(l  -  ?/m)(l  -  yh)  . 

Where  ?/£  is  the  total  cloud  fraction  over  a  grid  area. 
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Figure  2.4.  Definition  of  the  <7-layers  in  the  AFGL  GCM. 
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For  partly  cloudy  conditions,  the  radiative  heating/cooling  rate  at  each  model 
layer  is  obtained  by  linearly  weighting  the  percentages  of  the  total  cloud  cover,  r/(, 
and  clear  portion,  (1  —  ?/t),  in  the  form 

(f  r=*(§ 

where  the  superscripts  pc,  c,  and  nc  represent  partly  cloudy,  cloudy,  and  clear 
conditions,  respectively.  This  equation  applies  to  both  IR  and  solar  heating  rate 
calculations. 


29 


Section  3 


BULK  PARAMETERIZATIONS  OF  SOURCES  AND 
SINKS  OF  MICROPHYSICAL  PROCESSES 


Since  it  is  not  practical  in  a  large-scale  cloud  model  to  simulate  detailed  micro¬ 
physical  processes  for  each  individual  cloud  particle,  we  have  to  relate  these  pro¬ 
cesses  to  bulk  quantities,  which  represent  the  various  sources  and  sinks  for  cloud 
water  substances.  In  this  study,  the  parameterizations  of  these  bulk  quantities 
generally  follow  Kessler  (1969)  and  Sundqvist  (1978)  for  the  liquid  phase  and  Starr 
and  Cox  (1985),  Koenig  and  Murray  (1976),  and  Heymsfield  and  Platt  (1984)  for 
the  ice  phase. 

3.1  Liquid  Phase 

3.1.1  Condensation 


Condensation  transforms  water  vapor  into  cloud  droplets.  The  condensation 
occurs  as  a  result  of  saturated  expanding  air  and  may  be  related  to  the  time  rate  of 
change  of  the  saturation  vapor  pressure.  Consequently,  the  condensational  heating 
rate  may  be  derived  from  the  first  law  of  thermodynamics,  the  Clausius-Clapeyron 
equation  and  the  hydrostatic  approximation  as  follow: 


-  C,R,t 

■c  d  R,  ^CPRVT2  +  q,L2 


U)  . 


The  detailed  derivation  can  be  found  in  Appendix  B.  The  condensation  is  strongly 
dependent  on  the  vertical  velocity,  w.  As  saturated  air  risies,  i.e.,  w  >  0,  condensa¬ 
tion  takes  places  and  cloud  forms,  while  depletion  of  clouds  occurs  with  descending 
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saturated  air.  As  one  can  see  from  Equation  (2.1),  condensation  is  allowed  to  take 
place  at  a  relative  humidity  less  than  100%  but  greater  than  the  critical  relative 
humidity.  Otherwise,  cloud  cover  becomes  zero  and  no  condensation  occurs. 

3.1.2  Autoconversion 

Autoconversion  is  a  process  through  which  small  cloud  droplets  merge  into  large 
raindrops.  The  parameterization  of  autoconversion  follows  the  simple  exponential 
form  invented  by  Sundqvist  (1978),  as  follows  : 

P  =  C0qt{l  -  exp[-(qi/(Tjqir))2}}  ,  (3.1) 

where  Co  and  qir  are  parameters  denoting,  respectively,  the  typical  time  of  con¬ 
version  of  droplets  to  raindrops,  and  a  reference  value  of  cloud  LWC  above  which 
the  conversion  from  cloud  droplets  to  raindrops  increases  rapidly.  These  two  values 
act  a  way  similar  to  the  rate  of  conversion  and  the  parameterization  of  accretion 
discussed  in  Kessler  (1969). 

Kessler  proposed  two  separate  expressions  for  autoconversion  due  to  coalescence 
and  accretion  of  cloud  water  by  raindrops,  given  as  follows: 


.  =  Kx(qi  -  (Zi)  , 

(3.2) 

II 

tsT 

''•IS 

(3.3) 

where  qr  and  oq,  respectively,  are  the  rain  water  content  and  the  threshold  liquid 
water  content  below  which  there  is  too  little  cloud  water  to  produce  rain.  Figure  3.1 
demonstrates  a  comparison  of  the  autoconversion  rate  used  in  the  model  with  that 
proposed  by  Kessler(1969),  i.e.,  Equations  (3.2)  and  (3.3).  In  this  comparison,  we 
take  I<i  =  Co  =  10-3sec-I,ai  =  10~3,  K-i  =  2.19,  and  qr  =  5  x  10~5  (~  2mm/day). 
As  shown  in  Figure  3.1,  the  two  parameterizations  give  similar  results.  Apparently 
Sundqvist  (1978)  has  adequately  combined  the  above  two  equations  for  coalescence 
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Figure  3.1.  Comparison  of  the  present  parameterization  of  the  autoconversion 
rate  defined  in  Eq.  (3.1)  as  a  function  of  the  liquid  water  content  with  that  of 
Kessler(1969). 


and  accretion  into  a  simple  exponential  form.  This  form  has  been  successfully 
incorporated  into  the  GCMs  by  Sundqvist  (1981)  ,  and  Zheng  and  Liou  (1986). 

3.1.3  Evaporation  of  raindrops 

The  evaporation  of  raindrops  may  be  derived  from  the  combination  of  the  Fick’s 
law  of  the  mass  diffusion  and  the  Marshall-Palmer  raindrop  size  distribution.  The 
size  distribution  is  assumed  to  stay  the  same  throughout  the  whole  process  of 
evaporation.  The  detailed  derivation  of  raindrop  evaporation  is  shown  in  Appendix 
C,  and  the  formulation  is  given  by 

Er  =  I<e(l  -  h0)(  P  )042  ,  (3.4) 

where  the  evaporation  rate  constant  is 

I<e  =  5.95  x  10-3<k  /  . 

It  is  noted  that  the  ventilation  coefficient,  f,  is  included  in  the  calculation. 

Kessler  (1969)  followed  the  same  rules  but  used  a  different  approach  to  incor¬ 
porate  the  ventilation  effect.  He  derived  the  following  form  for  the  evaporation  of 
rain: 

Er  =  5.54  x  10 -%(pqr)°*5p-1  ,  (3.5) 

where  the  units  for  p  and  qi  are  in  gm/m3.  For  the  purpose  of  comparison,  the 
raindrop  evaporation  in  Equation  (3.5)  is  further  related  to  the  precipitation  flux 
as  follows: 

Er  =  K'e(l  —  h0)(  P  )0,55  ,  (3.6) 

Where  the  evaporation  rate  constant  is 

Kg  =  7.84  x  10"2?,  , 


Figure  3.2  shows  a  comparison  of  the  raindrop  evaporation  rates  in  Equations  (3.6) 
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PRECIPITATION  FLUX  (GM/CM8/SEC) 


Figure  3.2.  Comparison  of  the  present  parameterization  of  the  evaporation  rate, 
Er/q.l{l  —  h0 ),  based  on  Eq.  (3.4)  with  that  of  Kessler  (1969).  The  value  of  the 
ventilation  factor  used  in  the  present  parameterization  is  3. 


and  (3.4)  for  f=3.  The  differences  between  the  two  curves  are  quite  small  for  P{z) 
between  the  normal  range  of  large-scale  precipitation  of  10“2  and  10-1  gm/m2/sec. 

3.2  Ice  phase 

The  ice  phase  parameterizations  include  homogeneous/ heterogeneous  nucleation 
to  generate  ice  crystals,  depositional  growth  to  simulate  the  Bergeron-Findeison’s 
processes,  sublimation  and  melting  of  ice  crystals,  and  gravitational  settling  to 
deplete  the  ice  crystals.  The  Bergeron-Findeison’s  processes  occur  primarily  at 
temperatures  between  0°  and  —  40°  C,  while  homogeneous  nucleation  dominates 
below  —  40°  C.  The  parameterizations  of  these  source  and  sink  terms  in  the  ice 
phase  are  described  as  follows. 


3.2.1  Homogeneous  nucleation 


Homogeneous  nucleation  takes  place  whenever  the  condensation  occurs  at  tem¬ 
peratures  below  —  40°  C.  Homogeneous  nucleation  is  analogous  to  the  liquid  water 
condensation  process,  and  may  be  written  in  the  form 


Qi  =  Td 


<fj(T) .  LRa  —  CPRVT  .  _ 

~RT[CpRvP  +  qsP]W 


The  homogeneous  nucleation  rate  is  a  function  of  temperature  and  vertical  velocity. 


3.2.2  Heterogeneous  nucleation 

Heterogeneous  nucleation  forms  ice  crystals  on  active  ice  nuclei  primarily  at 
temperatures  between  0°  and  —  40°  C.  Based  on  the  experimental  work  of  Mason 
(1971),  Mason  and  van  der  Heuvel  (1959),  and  Koenig  and  Murray  (1976),  it  is 
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assumed  that  the  active  ice-forming  nuclei  (IN)  concentration  may  be  expressed  in 
terms  of  temperature  as  follow: 

N  =  A\exp{-^~-max[(T  —  To),  T*]}  , 

A2 

where  A\  =  1,  A2  =  4,  T*  =  0°  C,  and  T0  is  a  threshold  temperature,  which  is  taken 
to  be  em  —  20°  C.  Here,  it  is  assumed  that  when  the  saturated  air  temperature  is 
lower  than  the  T0,i.e.,  —  20°  C,  the  total  number  of  IN  remains  the  same  as  that 
at  T  =  —  20°  C,  as  shown  in  Figure  3.3. 

The  heterogeneous  nucleation  rate  of  ice  crystals,  Sh ,  may  be  given  by 

c  dN 

Sk  =  No¬ 


where  m0  is  the  mass  of  a  newly  nucleated  ice  crystal,  prescribed  as  10"11  gram, 
which  is  equivalent  to  an  ice  sphere  with  a  diameter  of  2.77/zm  (Koenig  and  Murray, 
1976),  and 


dN_ 

dt 


’  if  T  >  — 20°C 

,0  if  T  <  — 20°C 


Here,  Nt+&t  and  Nt  are  the  IN  concentrations  at  consecutive  time  steps  of 
and  t,  respectively. 


3.2.3  Bergeron-Findeison’s  processes 

The  Bergeron-Findeison’s  processes  transform  cloud  water  into  cloud  ice  at 
temperatures  between  0°  and  —  40°  C.  They  are  based  on  the  depositional  growth 
of  cloud  ice  at  the  expense  of  cloud  water.  The  depositional  rate  of  ice  crystals 
follows  that  developed  by  Koenig  (1971). 

For  a  single  ice  crystal,  Koenig  (1971)  derived  an  equation  for  the  depositional 
growth  of  ice  based  on  the  mass  diffusion  law  and  the  electrostatic  analogy  (Hallett, 
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1965).  The  growth  rate  is  a  function  of  temperature  because  of  the  strong  temper¬ 
ature  dependence  of  saturation  vapor  pressure  over  water  and  ice.  Koenig  showed 
that  the  deposition  rate  of  ice  crystals  as  a  function  of  temperature  has  two  maxima 
at  —  5°  and  — 15°  C.  This  result  agrees  qualitatively  with  experimental  values 
derived  by  Byers  (1965). 

To  simplify  the  complicated  vapor  diffusion  equation,  Koenig  (1971)  developed 
a  simple  equation  using  two  temperature-dependent  parameters  to  described  the 
crystal  vapor  diffusion  growth  rate  as  a  function  of  mass.  The  equation  is  written 
as  follows: 


dm 

dt 


=  a3m 


at 


(3.7) 


where  a3  and  a4  are  temperature-dependent  parameters  listed  in  Table  (3.1),  and 
m  is  the  mass  of  the  crystal  in  grams.  Figure  3.4  illustrates  the  curves  of  a3  and  a4 
as  functions  of  temperature.  The  maxima  at  —  5°  and  —  15°  C,  noted  previously, 
are  also  present  on  these  curves. 

The  bulk  quantity  of  the  ice  crystal  depositional  rate,  denoted  as  Sd,  may  be 
derived  from  Equation  (3.7)  by  multiplying  both  sides  by  N,  as  follows: 

(  Tt(rti)  =  3(  pq{  )04 

(3.8) 


{  fi  =  Qi/{qi  +  ft) 

where  /;  is  the  cloud  liquid  fraction  serving  as  an  adjustment  factor,  which  allows  for 
more  cloud  water  resulting  in  higher  ice  crystals  growth  rates.  In  addition,  growth 
is  stopped  if  there  is  no  cloud  water  to  supply  water  vapor.  Figure  3.5  illustrates  the 
ice  crystal  depositional  rate  as  a  function  of  temperature  for  different  values  of  cloud 
IWC(/5<7;).  It  is  noted  that  the  maximum  depositional  rate  occurs  at  —  15°  C,  and 
the  local  maximum  is  at  —  5°  C,  and  thelocal  minimum  is  at  —  8°  C.  As  a  matter 
of  fact,  the  maximum  growth-occurred  at  —  15°  C  is  due  to  a  large  vapor  pressure 
difference  and  a  large  capacitance  value  at  this  temperature.  This  maximum  ice 
growth  at  —15°  G  is  evident  in  our  numerical  results  shown  in  Section  5. 
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Table  3.1.  Values  of  parameters  a 3  and  a4  as  functions  of  temperature. 


TEMPERATURE(°C  ) 


Figure  3.5.  Ice  crystal  depositional  rate  as  a  function  of  temperature  for  three 
different  cloud  IWCs. 


3.2.4  Sublimation  and  melting  of  ice  crystals 

The  parameterization  of  the  rate  of  sublimation  is  based  on  the  work  proposed 
by  Koenig  and  Murray  (1976)  in  the  form 

S3  =  0.1037  (p  qa)  (1  -  h0)  ( p  g;)0,65  •  (3-9) 

This  expression  is  similar  to  the  expression  for  the  evaporation  of  rains  denoted 
in  Equation  (3.4).  Figure  3.6  demonstrates  the  ice  crystal  sublimation  rate  as  a 
function  of  cloud  IWC (pqi)  for  three  different  temperatures  with  ho  =  0.85. 

At  temperatures  above  0°C,  cloud  ice  crystals  are  assumed  to  melt  instanta¬ 
neously  back  to  cloud  water. 


3.2.5  Gravitational  settling 


Gravitational  settling  is  an  important  sink  for  cloud  IWC.  Therefore,  it  is  desir¬ 
able  to  have  a  representative  (inean-volume-weighted)  downward  flux  of  ice  water. 
Following  Starr  and  Cox  (1985),  the  downward  ice  flux  may  be  expressed  as 


nmv  dL  , 


(3.10) 


where  n,  m,  and  v  are  the  ice  crystal  size  distribution,  mass  and  terminal  velocity 
of  an  ice  crystal,  respectively,  and  L  is  the  maximum  dimension  of  an  ice  crystal. 
Lmin  and  Lmax  are  the  minimum  and  maximum  sizes  of  the  crystals  present.  The 
cloud  IWC  may  be  further  related  to  n  and  m  as 


P<ii 


nm  dL  . 


(3.11) 


In  order  to  calculate  the  downward  ice  flux,  the  size  distribution(n),  mass(m), 
and  terminal  velocity(u)  for  the  precipitating  ice  crystals  need  to  be  specifically 
expressed  in  terms  of  the  maximum  dimension  of  the  ice  crystals,  i.e.,  L. 

The  size  distributions  of  ice  crystals  in  mid-latitude  cirrus  clouds  were  measureed 

and  documented  by  Heymsfield  (1977)  using  an  optical  probe.  Based  on  those 
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Table  3.2.  Parameterized  coefficients  for  cirrus  particle  size  distribution. 
[After  Heymsfield  and  Platt  (1984)  ] 


Temperature  range 

an 

bn 

Mean  IWC 

-20°  ~  -25° 

uniT 

-2.56 

— iliWrgffflS— 

-25°  ~  -30" 

-2.51 

-30°  -  -35° 

mrrmm 

-2.21 

-35°  -  —40° 

WKT^rfmm 

-2.29 

-40°  ~  -45° 

-3.23 

-45°  ~  -50° 

mrrTiT— 

-3.15 

m^rrrrrmi 

-50°  ~  -55° 

-3.83 

WKEMMEE&M 

-55°  ~  -60° 

■m 

-3.85 

data,  Heymsfield  and  Platt  (1984)  further  averaged  all  spectra  of  crystals  within 
every  5°C  temperature  interval  ranging  from  —  20°C  to  —  6C°C.  Therefore,  the  ice 
crystal  size  distributions  in  terms  of  the  maximum  dimension  of  ice  crystals  were 
parameterized  for  different  temperature  ranges. 

The  size  distributions  are  written  in  the  form 

n  =  anLbn  ,  (3.12) 

where  n  is  the  number  concentration  (  )  and  L  is  the  particle  dimension 

( fim  ).  The  an  and  bn  are  the  intercept  parameter  and  the  slope,  respectively,  which 
are  temperature  dependent  and  given  in  Table  (3.2).  As  a  first  step,  the  parameters 
an  and  bn  are  obtained  from  the  temperature  range  of  —  30°C  ~  — 35°C,  since  the 
size  distributions  remained  rather  similar  in  the  temperature  range  —  20°  C  to 
-  40°  C. 

The  minimum  length  that  the  optical  probe  can  measure  is  about  20  jim.  Small 
ice  crystals  less  than  20  fim  sizes  could  be  missed  by  the  measurement  technique. 
Consequently,  small  ice  crystal  data  were  excluded  from  the  analysis  by  Heymsfield 
and  Platt  (1984).  For  this  reason,  20  //m  is  taken  as  the  minimum  length  of  ice 

crystals, jCmt-n,  which  is  needed  in  the  calculation  of  Equations  (3.10)  and  (  3.11  ). 
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Table  3.3.  Parameterized  coefficients  for  mass  of  an  ice  crystal. 


[After  Heymsfield  (1972)  ] 


Particle  habit 

bm 

Columns 

1.35  x  10"10 

1.7 

Bullet  rosettes 

4.4  x  10"14 

3.0 

Plates 

fffTTpia 

2.5 

3. 2. 5.1  Crystal  mass  and  terminal  velocity 

Heymsfield  (1972)  has  derived  empirical  expressions  relating  the  maximum  di¬ 
mension  of  an  ire  crystal  (L)  to  its  corresponding  mass  and  terminal  velocity.  Starr 
and  Cox  (1985)  successfully  incorporated  these  expressions  into  the  calculation 
of  downward  ice  flux  due  to  gravitational  settling  in  the  cirrus  cloud  model,  the 
purpose  of  which  was  partly  to  develop  cirrus  clouds  parameterizations  for  use  in 
large-scale  cloud  models.  In  this  study,  the  parameterizations  of  a  crystal  mass  and 
terminal  velocity  follow  those  used  in  Starr  and  Cox  (1985). 

The  individual  ice  crystal  mass,  m,  may  be  parameterized  in  terms  of  the 
maximum  dimension  of  an  ice  crystal,  L,  as  follows 

m  =  am  Lbm  ,  (3.13) 

where  m  and  L  are  in  units  of  grams  and  fim,  respectively.  The  constants  am  and 
bm  associated  with  various  ice  crystal  habits  are  given  in  Table  (3.3). 

The  terminal  velocity  v  in  ms-1  of  an  ice  crystal  of  maximum  dimension  L  in 
lim  is  given  by 

v  =  avLb"  .  (3.14) 

The  constants  av  and  bv  are  dependent  on  the  ice  crystal  habit  and  size  as  given  in 

Table  (3.4).  These  coefficients  were  derived  at  an  ambient  pressure  of  400  mb.  For 
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Table  3.4.  Parameterized  coefficients  for  terminal  velocity  of  an  ice  crystal 


Size  range 
(/im) 


0-200 


200-400 


400-600 


800-5000 


Size  range 
(urn) 


0-200 


200-400 


400-600 


600-800 


800-5000 


Size  range 
(/zm) 


200-400 


400-600 


600-800 


800-5000 


[After  Starr  and  Cox  (1985)  ] 


Columns 


an  ice  crystal  at  a  given  pressure  p,  Beard  and  Pruppacher  (1969)  suggested  that 
its  terminal  velocity,  u(p),  may  be  given  by 


«(p)  =  «(?*)  (p*/p)’ 

where  p*  denotes  a  reference  pressure  level,  and  v{p*)  is  the  terminal  velocity  of  a 
crystal  at  the  reference  level.  To  obtain  the  terminal  velocity  at  any  given  pressure, 
we  simply  apply  the  pressure  adjustment  to  that  derived  from  Equation  (3.14). 

3. 2. 5. 2  Downward  ice  flux 

Since  the  parameterizations  of  the  ice  crystals  size  distribution  (n),  ice  crystal 
mass  (m),  and  terminal  velocity  (v)  are  all  parameterized  in  terms  of  the  maximum 
dimension  (L),  the  downward  ice  flux  may  be  obtained  by  analytically  integrating 
Equation  (3.10),  provided  that  the  size  of  the  largest  crystals,  Lmax,  is  known. 
Since  the  lower  limit,  Lm,-n,  is  set  to  be  20  pm ,  as  stated  previously,  the  upper  limit 
{Lmax)  may  be  determined  from  Equation  (3.11)  in  the  following  way.  Once  pqi 
for  the  n-th  time  step  is  obtained,  Equations  (3.12)  and  (3.13)  are  substituted2 into 
Equation  (3.11)  to  calculate  Lmax.  Having  determined  Lmax ,  the  downward  ice  flux 
can  be  computed  from  Equation  (3.10)  as  follows: 

Mp*)  =  .  . Yrr.  - iL+‘"+‘,+1)]  (3.i5) 

The  downward  flux  is  the  function  of  Lmax,  which  in  turn  is  the  function  of  cloud 
IWC  {pqi).  Consequently,  the  downward  ice  flux  may  be  expressed  in  terms  of  cloud 

IWC. 

Figure  3.7  demonstrates  the  computed  ice  downward  flux  as  a  function  of  cloud 
IWC  for  the  temperature  range  of  —30°  ~  —35°  C.  In  this  calculation,  all  crystals 
were  assumed  to  be  columns.  When  cloud  IWCs  varies  from  5  x  10-2gm/m3, 
the  ice  precipitation  flux  increases  from  2  x  10-4  to  about  10“2gm/m2/sec,  which 
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corresponds  to  about  6.3  to  310  mm/yr.  Comparing  these  values  of  ice  precipitation 
with  the  annual  global  precipitation  of  1000  mm  estimated  by  Jaeger  (1983),  we 
note  that  gravitational  settling  is  a  significant  process  in  the  depletion  of  the  IWC 
and  the  production  of  ice  precipitation. 
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Section  4 


MODEL  INPUT  AND  VERIFICATION  DATASETS 
4.1  Model  input  datasets 

The  analysis  atmospheric  data  generated  from  the  AFGL  spectral  global  model 
at  the  12Z  of  July  1-4,  1979,  are  used  as  input  in  the  cloud  model.  These  data 
contain  temperature,  specific  humidity,  and  wind  fields  (t2,u,  u))  at  all  12  a  — layers 
in  the  atmosphere  plus  the  surface  pressure  and  geopotential  height.  These  analysis 
data  were  carried  out  at  a  rhombodial  15  (R15)  truncation  which  corresponds  to  a 
horizontal  resolution  of  38  and  48  grid  points  in  the  latitudinal  and  longitudinal, 
respectively. 

In  the  vertical,  since  these  GCM-generated  input  data  are  available  at  the  a  — 
coordinate  shown  in  Figure  2.4,  it  is  necessary  to  perform  vertical  interpolations  to 
obtain  these  values  in  the  z-coordinate.  A  logarithmic  interpolation  was  performed 
on  the  AFGL  atmospheric  profiles  to  map  them  from  the  a  levels  to  the  height  levels 
of  the  cloud  model.  All  variables  except  temperature  were  converted  in  this  manner. 
The  temperatures  were  converted  using  a  linear  interpolation.  Furthermore,  the 
vertical  velocity  generated  from  the  AFGL  spectral  model  is  in  a  a  —  coordinate 
and  is  denoted  as  a.  The  conversion  of  cr  to  w  is  accomplished  based  on  the  following 
relationship  suggested  by  Kasahara(1974). 

,dz .  _  _  .  dz 

u.  =  (^  +  vff.v„,+^ 

where  z  is  the  geopotential  height. 
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Figures  4.1  -  4.5  show  the  zonal  mean  values  of  temperature,  specific  humid¬ 
ity,  and  wind  fields  (u,  u,  w )  averaged  over  July  1-4,  1979.  As  one  can  see  from 
Figure  4.1,  the  maximum  tropospheric  temperature  occurs  near  10°  N.  The  tro¬ 
pospheric  temperatures  within  the  Antarctic  Circle  are  much  lower  than  those 
within  the  Arctic  Circle.  The  mean  temperature  distribution  in  both  hemispheres 
is  characterized  by  an  isobarically  horizontal  decrease  poleward  at  all  levels  of  the 
troposphere  except  above  about  150  mb  in  the  tropics,  where  the  tropopause  is 
the  highest.  Due  to  the  lack  of  solar  radiation  within  the  Antarctic  Circle  in  July, 
the  overall  pole-to-equator  temperature  contrast  in  the  troposphere  is  much  greater 
in  the  Southern  Hemisphere  than  in  the  Northern  Hemisphere.  It  is  noted  that  a 
temperature  inversion  occurs  in  the  Northern  Hemisphere  polar  region. 

Figure  4.5  shows  the  zonal  average  of  vertical  velocity  in  the  z-coordinate,  i.e.,  w. 
The  rising  motions  are  associated  with  the  equator  and  the  middle  latitudes  in  both 
hemispheres.  In  the  vicinity  of  the  equator,  the  upward  motions  are  stronger  in 
the  Northern  Hemisphere  than  in  the  Southern  Hemisphere.  The  strongest  vertical 
velocity  occurs  at  the  latitude  near  10°  N  in  the  branch  of  rising  motion  of  the 
Northern  Hemisphere  Hadley  cell.  The  downward  motions  are  found  in  the  vicinities 
of  the  subtropic  highs,  i.e.,  30°  N  and  20°  S.  The  downward  motions  in  the  Southern 
Hemisphere  are  stronger  than  those  in  the  Northern  Hemisphere. 

The  radiative  transfer  model  also  requires  values  for  the  surface  albedo,  which 
is  defined  as  the  portion  of  the  incoming  solar  radiation  that  is  reflected  by  the 
earth’s  surface.  A  global  dataset  of  surface  albedo  for  July  1979  has  been  compiled 
by  Griffin  (1987)  and  Mitchell  (1990).  This  dataset,  shown  in  Figure  4.6,  is  used  in 
this  study.  As  one  can  see,  the  surface  albedo  can  vary  in  value  from  a  few  percent 
(ocean  surface)  to  more  than  80%  (snow  or  ice  cover)  with  sharp  changes  occurring 
along  the  continental  coastlines.  The  surface  albedo  is  larger  in  the  antarctic  area 
than  in  the  arctic  region. 
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Figure  4.3.  The  analysis  field  of  zonal  wind  averaged  over  July  1-4,  1979 


Figure  4.4.  The  analysis  field  of  meridional  wind  averaged  over  July  1-4,  1979. 
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Figure  4.6.  Monthly  averaged  surface  albedo  for  July,  1979. 


4.2  Verification  datasets 

It  is  one  of  the  main  objectives  in  this  study  to  assess  the  model  performance 
by  verifying  the  predicted  results  against  the  verification  datasets.  Ideally,  such 
datasets  should  provide  simultaneous  measurements  of  cloud  cover,  cloud  LWC 
and  IWC,  and  precipitation  rates  on  a  global  basis. 

Although  ground-based  cloud  cover  observations  have  been  recorded  for  decades, 
the  cloud  amount  obtained  by  a  ground  observer  is  subjective,  tends  to  be  overesti¬ 
mated  (Hughes,  1984)  and  is  less  reliable  at  night.  For  multilayer  cloud  conditions, 
the  cloud  amounts  associated  with  the  upper  cloud  layers  are  frequently  under¬ 
estimated  due  to  the  lower  clouds  obscuring  the  upper  clouds.  It  is  the  greatest 
intrinsic  weakness  of  the  ground-based  cloud  observations  that  most  observational 
stations  are  unevenly  gathered  in  high  population  areas  and  can  provide  only  sparse 
observational  data  compared  to  satellite  observations. 

In  contrast  to  surface-based  observations,  satellite  observations  are  objective 
and  the  infrared  sensors  on  board  satellites  can  provide  nightime  observations 
of  as  good  a  quality  as  those  taken  in  daytime.  Satellite  observations  tend  to 
underestimate  the  cloud  amounts  associated  with  the  lower  cloud  layers,  which 
are  obscured  by  the  higher  cloud  layers.  In  addition,  there  is  substantial  error 
arising  from  the  failure  to  distinguish  between  snow/ice  surfaces  and  cloud  tops, 
especially  in  the  polar  regions.  Nevertheless,  satellite  observations  provide  an  even 
global  scope  of  cloud  cover  with  resolution  as  high  as  48  km  by  48  km  per  pixel. 
This  high  resolution  global  coverage  of  satellite  observed  cloud  cover  makes  satellite 
observations  superior  to  other  cloud  observations. 

4.2.1  3DNEPH  cloud  data 

The  Air  Force  Global  Weather  Central  (AFGWC)  high  res  mtion  three  dimen¬ 

sional  ncphanalysis  program,  known  as  3DNEPH,  is  the  only  operational  cloud 
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analysis  that  combines  both  satellite  and  conventional  reports.  The  3DNEPH 
cloud  information  is  derived  from  the  DMSP  (Defense  Meteorological  Satellite 
Program)  satellites,  which  are  sun-synchronous  polar  orbiting  satellites  with  visual 
and  infrared  sensors. 

3DNEPH  generates  high  resolution  global  cloud  information  every  3  and  6  hours 
for  the  entire  Northern  and  Southern  Hemispheres,  respectively.  Each  hemisphere  is 
divided  into  64  3DNEPH  boxes  as  illustrated  in  Figure  4.7,  with  each  box  containing 
4096  (64x64)  grid  points.  The  horizontal  resolution  of  3DNEPH  is  approximately 
48  km  at  60°  latitude.  The  vertical  grid  consists  of  15  layers  of  variable  thickness 
which  extend  from  the  surface  up  to  17  km.  3DNEPH  cloud  analysis  consists  of 
cloud  amount,  cloud  base,  cloud  top,  and  low,  middle,  high  cloud  types. 

Based  on  3DNEPH  analysis,  the  global  cloud  climatologies  for  January  and  July 
1979  were  processed  into  appropriate  formats  for  model  verification.  This  work 
had  been  done  mostly  by  Koenig  (1985)  and  Griffin  (1987),  and  partly  by  Mitchell 
(1990).  The  global  distribution  of  cloud  amount  on  July  1-4,  1979  is  derived  from 
the  above  cloud  climatologies  and  used  as  the  verification  dataset  for  the  cloud 
cover  simulated  in  the  model.  In  this  dataset,  the  original  3DNEPH  15-layer  cloud 
amount  is  combined  to  produce  low,  middle,  high,  and  total  (LMHT)  cloud  data 
bases.  The  low/middle  cloud  boundary  is  defined  as  terrain  height  -1-  2  km,  while 
the  middle/high  cloud  boundary  is  terrain  height  +  the  height  of  the  —  20°  C 
isotherm.  In  addition,  the  high  resolution  (48  km)  3DNEPH  horizontal  grid  points 
were  mapped  to  the  coarser  Nimbus  7  ERB  sub-target  area  grid  structure  (160  km). 

It  is  important  to  recognize  the  limitation  of  3DNEPII  cloud  analysis.  A  thresh¬ 
old  method  is  employed  to  process  satellite  images  to  obtain  the  3DNEPH  cloud  in¬ 
formation.  The  visual  image  is  compared  with  background  brightness;  The  infrared 
image  is  compared  with  a  threshold  temperature,  which  is  largely  determined  by 
the  surface  temperature.  If  an  observed  brightness  temperature  of  a  pixel  is  colder 
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Figure  4.7.  Configuration  of  the  3DNEPII  grid  over  the  Northern  Hemisphere. 


than  the  threshold  temperature,  then  the  cloudy  condition  is  assigned  to  the  pixel. 
Consequently,  the  cloud  amount  tends  to  be  overestimated  when  the  associated 
surface  temperature  is  low,  and  underestimated  when  a  temperature  inversion  is 
located  in  the  region.  This  may  be  evident  from  the  fact  that  3DNEPH  over-  and 
underestimates  cloud  amounts  in  the  wintertime  and  summertime  polar  regions, 
respectively  (Henderson-Seller,  1986). 

The  high  cloud  amount  is  underestimated  in  both  infrared  and  visual  processes 
in  3DNEPH.  Satellite  infrared  sensors  tend  to  interpret  high  thin  nonblack  clouds  as 
middle  clouds  and  underestimates  the  high  clouds  amount.  Since  the  transmissivity 
of  these  clouds  is  nonzero,  infrared  energy  from  below  the  cloud  is  transmitted  by 
the  cloud  making  it  appears  warmer  than  it  actually  is.  The  visual  sensor  also 
fails  to  detect  small  cloud  amounts  in  the  vicinity  of  areas  of  strong  sun  glint. 
Frequently,  thin  cirrus  over  a  bright  background  (for  example,  over  a  desert)  is  not 
detected  by  the  visual  sensor. 

The  aforementioned  problems  become  evident  in  attempts  to  use  the  3DNEPH 
cloud  data  to  verify  the  predicted  cloud  cover  in  Section  5.  However,  several 
intercomparisons  of  the  cloud  cover  between  3DNEPH  and  other  cloud  retrieval 
techniques  have  been  made  (Koenig  and  Liou,  1983;  Gordon  etal.,  1984;  Hughes, 
1984),  and  they  indicated  that  the  3DNEPH  total  cloud  cover  appears  to  be  reliable, 
except  in  the  polar  regions  where  all  satellite  retrieval  methodologies  tend  to  fail. 
Alternatively,  the  surface  observational  cloud  climatology  compiled  by  London 
(1957)  is  adopted  in  this  study  as  a  complementary  verification  data  set  for  the 
polar  regions. 


4.2.2  Earth  Radiation  Budget  (ERB)  dataset 

The  Earth  Radiation  Budget  (ERB)  instrument  aboard  the  Nimbus  7  satellite 
takes  broadband  measurements  of  the  solar  and  terrestrial  radiation  reflected  and 
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emitted  by  the  earth-atmosphere  system  (Jacobowitz  etal.,  1978).  The  reflected 
solar  radiation  is  referred  as  albedo,  and  the  emitted  infrared  radiation  is  called 
outgoing  longwave  radiation  (OLR).  These  OLR  and  albedo  data  have  been  exam¬ 
ined  and  published  in  a  number  of  articles  (Jacobowitz  etal.,  1984  a,b;  Arking  and 
Vemury,  1984;  Kyle  and  Vasanth,  1986),  and  appear  to  be  an  ideal  earth  radiation 
data  base  for  use  in  the  verification  of  large-scale  cloud  models.  In  this  study,  the 
observed  albedo  is  converted  to  the  net  solar  flux  at  TOA  to  be  compared  with  the 
model-calculated  net  solar  flux. 

The  archived  ERB  flux  data  tapes  for  the  July,  1979  provided  by  the  NASA  ERB 
team  was  used  in  this  study  to  validate  the  earth  radiation  budget,  i.e.,  OLR  and 
albedo,  produced  by  the  model.  Since  they  are  strongly  modulated  by  the  radiative 
effects  of  clouds,  including  cloud  cover  and  cloud  LWC/IWC,  the  ERB  data  have 
been  widely  used  in  large-scale  cloud  model  verifications.  (Slingo,  1986;  Smith, 
1990).  In  particular,  when  utilized  together  with  3DNEPH,  the  ERB  data  may 
provide  valuable  guidance  on  the  cloud  radiative  properties  of  cloud  LWC/IWC 
simulated  in  the  model. 

Nimbus  7  is  a  sun-synchronous  polar  orbiting  satellite  with  equator  crossings  at 
local  noon  (ascending)  and  midnight  (descending).  Terrestrial  radiation  is  measured 
by  wide  field-of-view  (WFOV)  and  narrow  field-of-view  (NFOV)  radiometers.  The 
WFOV  channels  integrate  the  outgoing  radiation  over  the  entire  earth’s  disc  visible 
from  the  satellite;  the  NFOV  channels  measure  the  outgoing  radiance  from  a 
number  of  directions  with  a  mean  resolution  of  150  km  by  150  km. 

The  NFOV  radiance  observations  from  Nimbus  7  are  sorted  into  2070  Target 
Areas  (TAs),  which  cover  the  whole  globe  and  are  roughly  of  equal  area.  Each 
TA  is  approximately  500  km  x  500  km  and  is  bounded  by  latitude  and  longitude 
lines.  There  is  a  4.5°  interval  between  the  latitude  lines.  For  each  latitude  band, 
the  longitudinal  increments  for  a  TA  vary  from  4.5°  at  the  Equator  to  120°  near 
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the  pole.  Each  TA  is  further  divided  into  nine  (3x3)  Subtarget  Areas  (STAs).  Each 
STA  roughly  corresponds  to  an  area  of  160x160  km2.  The  exact  boundaries  for 
the  2070  TAs  are  jhown  in  Table  (4.1)  and  displayed  on  a  Mercator  projection  in 
Figure  4.8. 

To  derive  mean  daily  flux  values  of  the  earth  radiation  budget,  angular  depen¬ 
dence  models  and  diurnal  correction  models  are  adopted.  An  angular  dependence 
model  converts  the  NFOV  radiance  measurement  to  a  flux  at  TOA;  the  diurnal 
models  are  used  to  estimate  the  daily  averaged  albedo  from  a  once  a  day  observation. 
Errors  involved  in  the  angular  model  conversion  are  smaller  than  those  in  the  diurnal 
models.  Since  the  reflectivity  is  a  function  of  the  angle  of  incident  radiation,  the 
instant  albedo  measurement  is  not  a  good  estimate  of  the  mean  daily  albedo.  Such 
errors  were  carefully  examined  and  reduced  in  the  NASA  ERB  data  by  Kyle  etal., 
(1990). 


4.2.3  Scanning  Multichannel  Microwave  Radiometer  (SMMR) 

Monthly  mean  distributions  of  liquid  water  were  derived  from  the  Nimbus  7 
SMMR  observations  over  the  oceans  for  the  period  November  1978  to  November 
1979  (Prabhakara  and  Short  1984).  The  liquid  water  estimate  for  the  July,  1979 
is  applicable  to  the  monthly  average  over  an  area  approximately  400  km  on  a 
side  (3°  x  5°,  latitude/longitude).  The  liquid  water  derived  includes  vertically 
integrated  amounts  of  liquid  droplets  in  both  clouds  and  rain.  Since  the  ice-free  sea 
surface  emissivity  is  nearly  constant,  the  liquid  water  verification  data  were  limited 
to  the  oceans  equatorward  of  60°  latitude.  Over  land,  there  is  no  climatological 
liquid  water  data  on  the  global  extent  for  comparison. 

Liquid  droplets  in  clouds  and  rain  absorb  and  scatter  microwave  radiation. 
The  microwave  brightness  temperatures,  observed  by  the  SMMR,  are  significantly 
modulated  by  liquid  droplets  in  the  atmosphere.  These  droplets  produce  a  bright- 
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Table  4.1.  Nimbus-7  ERB  TA  boundaries  in  longitudinal  and  latitudinal  lines. 


Sequential 

Longitude 

Target  No. 

Latitude 

Limits 

Interval* 

Southern 

Northern 

Lower 

Upper 

Hemisphere 

Hemisphere 

Limit 

Limit 

956-1035 

1036-1115 

Equator 

4.5 

4.5 

876-955 

1116-1195 

4.5 

9.0 

4.5 

796-875 

1196-1275 

9.0 

13.5 

4.5 

716-795 

1276-1355 

13.5 

18.0 

4.5 

644-715 

1356-11)27 

18.0 

•  22.5 

5.0 

572-6143 

1 D28-1 D99 

22.5 

27.0 

5.0 

500-571 

1500-1571 

27.0 

31.5 

5.0 

1)28-1499 

1572-1643 

31.5 

36.0 

5.0 

368-1)27 

1644-1703 

36.0 

40.5 

6.0 

308-367 

1704-1763 

40.5 

45.0 

6.0 

21)8-307 

1764-1823 

45.0 

49.5 

6.0 

200-2D7 

1824-1871 

49.5 

54.0 

7.5 

155-199 

1872-1916 

54.0 

58.5 

8.0 

115-15D 

1917-1956 

58.5 

63.0 

9.0 

79-111) 

1957-1992 

63.0 

67.5 

10.0 

1)9-78 

1993-2022 

67.5 

72.0 

12.0 

29-1)8 

2023-2042 

72.0 

76.5 

18.0 

13-28 

2043-2058 

76.5 

81.0 

22.5 

1)— 1 2 

2059-2067 

81 .0 

85.5 

40.0 

1-3 

2068-2070 

85.5 

Pole 

120.6 

#For  each  latitude  band,  the  longitude  intervals  start  at  the  0° 


meridian  and  progress  west  by  the  increments  listed. 
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Figure  4.8.  Nimbus-7  ERB  global  grid  of  2070  target  areas. 


ness  temperature  change  that  increases  significantly  with  increasing  frequency. 
However,  the  Rayleigh  approximation  for  rain  drops  breaks  down  as  the  frequency 
increases  to  some  extent.  As  a  result,  SMMR  brightness  measurements  at  the  high 
frequency  channels,  18,  21  and  37  GHz,  do  not  increase  linearly  with  the  LWC  in 
the  atmosphere.  This  can  adversely  affect  the  remote  sensing  of  LWC.  For  this 
reason,  only  the  lower  frequency  channels,  6.6  and  10.7  GHz,  are  used  to  derive  the 
liquid  content  in  the  atmosphere. 

Since  there  is  no  ground  truth  with  which  to  compare,  the  accuracy  of  liquid 
water  observations  from  SMMR  is  roughly  estimated  as  ~  10  mg/cm2  based  on 
theoretical  considerations  (Prabhakara  etal.,  1983).  However,  by  taking  monthly 
averages,  this  accuracy  may  be  improved  to  5  mg/cm2  (Prabhakara  and  Short, 
1984).  The  temporal  and  spatial  variation  of  liquid  water  is  very  great.  For 
example,  while  instantaneous  LWCs  in  deep  precipitating  clouds  can  be  as  high  as 
3000  mg/cm2,  monthly  average  amounts  over  large  areas  rarely  exceed  30  mg/cm2. 
This  is  a  consequence  of  the  relatively  small  spatial  scales  and  short  lifetimes  of 
cloudiness.  In  addition,  the  SMMR  observations  seem  to  underestimate  cloud  LWC. 
Cloud  IWC  is  transparent  to  microwaves  at  the  SMMR  frequencies. 

The  model  input  and  verification  datasets,  including  AFGL  analysis  data,  ERB, 
3DNEPH  and  SMMR,  are  summerized  in  Table  (4.2). 


Table  4.2.  A  summerized  table  for  model  input  and  verification  datasets. 
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Section  5 


MODEL  PERFORMANCE  AND  VERIFICATION 
5.1  Design  of  the  Verification 

In  the  development  of  any  parameterization  scheme,  the  prediction  results  must 
be  carefully  verified  to  assess  the  performance  of  the  new  scheme.  A  schematic 
illustration  of  the  model  verification  is  given  in  Figure  5.1.  The  large-scale  cloud 
model  was  integrated  for  96  hours  with  a  time  step  of  30  mintues  from  the  initial 
conditions  taken  at  12Z  of  July  1,  1979.  The  winds  to  circulate  cloud-moisture 
are  not  predicted  in  the  cloud  model.  Instead,  they  are  prescribed  as  four  sets  of 
analysis  wind  fields  on  each  12Z  of  July  1  to  4,  1979.  During  the  96-hour  model 
integration,  each  set  of  12Z  wind  fields  is  employed  at  each  time  step  for  24  hours. 
It  is  noted  there  is  no  cloud-radiation  feedback  to  the  dynamic  structure  in  this 
study. 

With  input  of  the  initial  data  and  wind  fields,  the  cloud  model  is  capable  of  fore¬ 
casting  large-scale  cloud  cover  and  cloud  LWC/IWC,  which  are  interactively  input 
into  the  radiation  model.  The  radiation  model  calculates  radiative  heating/cooling 
rates  based  on  the  predicted  cloud  properties.  The  radiative  heating/cooling  rates, 
in  turn,  feed  back  to  the  cloud  model  through  the  thermodynamic  equation,  and 
ultimately  modify  the  cloud  fields.  Therefore,  cloud-radiation  interaction  through 
the  medium  of  the  cloud  cover  and  cloud  LWC  and  IWC  is  accomplished  with  the 
aid  of  the  cloud  model.  In  addition,  the  radiation  model  computes  OLR  and  net 
solar  flux  at  TOA  to  be  compared  with  the  ERB  satellite  observational  data.  In 
this  study,  the  interactive  radiation  computations  are  undertaken  every  three  hours, 
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Figure  5.1.  A  schematic  illustration  of  the  model  verification. 


while  the  calculated  heating/cooling  rates  are  input  in  the  cloud  model  at  every 
time  step  during  those  three  hours. 

The  simulated  cloud  cover,  earth’s  radiation  budget  and  cloud  LWC  and  IWC 
are  verified  against  3DNEPH  cloud  cover,  ERB  and  SMMR,  as  well  as  in  situ  ob¬ 
servations  and  other  model; studies.  However,  there  are  some  problems  which  make 
the  use  of  satellite  observations  in  model  verification  not  entirely  straightforward. 
For  instance,  verification  datasets  such  as  3DNEPH  and  ERB  are  based  on  the 
daily  average  over  local  time,  while  the  model  simulation  is  evolved  in  Greenwich 
Mean  Time  (GMT).  Allowing  intercomparisons  to  be  made  between  the  satellite 
observations  and  model  simulations,  we  have  to  average  the  model  predicted  fields 
over  local  time  in  order  to  obtain  daily  averaged  model  results.  Unless  specified, 
all  of  the  results  including  satellite  observations  and  model  simulation  are  averaged 
over  July  3-4,  1979.  That  is  most  of  the  results  shown  in  this  and  latter  section  are 
presented  as  a  2-day  average  of  July  3  and  4. 


5.2  Initial  Data 

The  analysis  data  on  12Z,  July  1,  1979  for  temperature,  specific  humidity,  and 
winds (ft,u,u>),  shown  in  Figures  5.2  -  5.6,  are  used  as  initial  values  to  perform 
a  96-hour  model  simulation.  Data  for  the  cloud  LWC/IWC  are  not  available  for 
initial  values  and  LWC/IWG  are  set  as  zero  initially. 

The  initial  cloud  cover  may  be  deduced  from  the  initial  temperature  and  specific 
humidity.  The  initial  temperature  is  used  to  derive  the  saturation  specific  humidity, 
denoted  by  qs(T),  which  is  defined  as  follows: 

h,(t )  =  e-^p- 

where  e  =  0.622  cu.d  es(T)  is  a  saturation  vapor  pressure,  which  is  related  to 
large-scale  temperature  by  the  Clausius-Clapeyron  equation 

e,(T)  =  0.622  X  exp(l  -  -L) 
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Figure  5.3.  The  zonal 
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Figure  5.4.  The  zonal  mean  field  of  initial  zonal  wind. 
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Figure  5.6.  The  ^onai  mean  field  of  initial  vertical  velocity. 


Since  saturation  vapor  pressure  is  a  function  of  temperature  only,  saturation  specific 
humidity  is  determined  by  temperature  and  pressure.  Initial  relative  humidity,  the 
ratio  of  initial  specific  humidity  to  its  saturation  value,  is  determined.  Therefore, 
the  initial  cloud  cover  is  obtained  through  the  definition  in  Equation  (2.6  )  with 
the  values  of  h0  derived  from  the  Smagorinsky’s  (1960)  empirical  expression  shown 
in  Fig.  2.1. 

The  zonal  average  of  initial  cloud  cover  is  shown  in  Figure  5.7.  It  is  noted  that 
clouds  may  be  formed  in  every  model  layer,  except  the  lowest  three  and  the  top 
two  layers.  One  of  the  distinctive  features  in  the  figure  is  the  small  cloudiness 
associated  with  the  subtropic  highs  at  30°  N  and  20°  S,  where  downward  motions 
are  dominant.  Large  cloudiness  is  found  at  the  tropics,  the  middle  latitudes  in 
the  both  hemispheres,  and  the  arctic  regions.  The  large  cloudiness  is  associated 
with  upward  motions  in  the  intertropical  convergence  zone  (ITCZ)  and  mid-latitude 
storm  tracks.  The  fraction  of  high  cloud  is  about  0.2;  low  and  middle  cloud  covers 
are  larger  than  high  cloud  cover.  The  cloudiness  within  the  Antarctic  Circle  is  less 
than  that  within  the  Arctic  Circle. 

5.3  Model  Results  and  Verification 

5.3.1  Zonal  means 

5. 3. 1.1  Zonal  mean  cloud  cover 

In  Figure  5.8,  the  model-predicted  zonal  mean  total  cloudiness,  averaged  over 
July  3-4,  is  compared  with  the  corresponding  3DNEPH  cloud  observations.  Dif¬ 
ferences  between  the  two  curves  are  generally  within  10%  of  the  total  cloud  cover 
except  in  the  tropics,  arctic,  and  antarctic  areas. 

As  was  mentioned  in  the  Section  4,  3DNEPII  under-  and  over-estimates  cloud 
covers  in  arctic  and  antarctic  regions,  respectively,  for  the  summertime  of  the 
Northern  Hemisphere.  This  observational  bias  is  depicted  in  Figure  5.9  by  the 
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Figure  5.8.  Comparisons  of  the  zonal  mean  total  cloudiness  in  the  simulation  and 
3DNEPH. 
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comparison  of  3DNEPH  monthly  mean  data  with  London’s  cloud  climatology  based 
on  surface  cloud  observations.  Since  London’s  cloud  climatology  was  compiled 
only  for  the  Northern  Hemisphere,  the  comparison  is  confined  to  the  Northern 
Hemisphere.  As  indicated  in  Figure  5.9,  London’s  cloud  climatology  shows  about 
70%  of  total  cloudiness  in  the  summertime  Arctic  region  and  40%  in  the  wintertime. 
The  70%  of  total  cloudiness  in  the  Arctic  is  close  to  our  model  prediction  in  the 
area.  This  explains  the  differences  between  the  predicted  and  3DNEPH  cloud  cover 
in  the  arctic  and  antarctic  regions. 

The  observed  large  cloudiness  in  the  tropics  is  largely  contributed  by  tropical 
anvils,  which  originate  from  convective  clouds.  Since  the  present  large-scale  cloud 
model  is  developed  primarily  for  stratiform  clouds,  the  simulated  cloud  cover  in  the 
tropics  is  underestimated,  as  is  evident  in  Figure  5.8.  Another  shortcoming  in  the 
simulated  cloud  cover  is  the  overestimation  of  cloud  cover  near  the  60°  S,  where 
the  incoming  solar  radiation  is  small.  As  a  result,  longwave  radiation  produces 
strong  radiative  cooling  at  the  top  of  middle  clouds  and  results  in  a  large  increase 
of  total  cloudiness  in  the  region.  It  appears  that  improvement  could  be  made  if 
a  dynamic  instability  feedback  were  allowed  in  the  simulation.  Nevertheless,  the 
comparison  of  zonal  mean  cloud  cover  shows  good  agreement  between  the  simulation 
and  observation.  The  simulated  small  cloud  cover  located  at  the  subtropic  highs 
in  both  of  the  hemispheres  are  realistic.  A  large  cloudiness  in  the  tropics  is  also 
predicted,  although  the  3DNEPH  data  suggests  that  it  is  underestimated. 

5. 3. 1.2  Zonal  mean  earth’s  radiation  budget 

Figures  5.10  and  5.11  are  the  same  as  Figure  5.8  except  for  OLR  and  net  solar 
flux  at  TOA,  respectively.  Figure  5.10  indicates  good  agreement  in  OLR  between 
the  model  simulation  and  ERB  observation  with  differences  being  generally  less 
than  10  W/m2.  Both  curves  indicate  two  maxima  located  in  the  subtropics  in 
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Figure  5.9.  Comparisons  of  the  3DNEPH  and  London’s  cloud  climatology. 


both  hemispheres  with  a  minimum  in  the  tropics.  These  two  well-defined  OLR 
maxima  are  associated  with  the  subtropical  highs,  where  small  cloudiness  and  high 
temperature  produce  strong  outgoing  longwave  flux.  The  OLR  in  the  Southern 
Hemisphere  is  considerably  smaller  than  in  the  Northern  Hemisphere  because  of 
the  low  temperatures  in  the  wintertime  Southern  Hemisphere. 

The  minimum  OLR  in  the  tropics  is  related  to  the  ITCZ  region,  in  which  anvil 
cirrus  originating  from  deep  convective  clouds  occupy  a  large  part  of  the  sky  at 
very  high  altitudes.  These  anvil  clouds  absorb  the  OLR  emitted  from  the  earth’s 
surface  and  re-emit  to  space  at  much  lower  temperature  than  at  the  surface.  As 
a  consequence,  the  OLR  in  the  ITCZ  is  greatly  reduced.  The  minimum  OLR  in 
the  tropics  is  reasonably  well-predicted,  although  the  simulated  value  is  30  W/m2 
higher  than  that  in  the  ERB  data.  This  overestimation  of  OLR  may  be  explained  by 
the  following  two  reasons.  First,  the  altitude  of  high  clouds  defined  in  the  model  is 
usually  lower  than  the  actual  height  of  a  tropical  anvil  cloud.  Second,  the  simulated 
cloud  covers  in  the  tropics  are  underestimated. 

Illustrated  in  Figure  5.11  is  the  zonal  mean  net  solar  downward  flux  at  TO  A 
for  the  model  simulation  and  the  ERB  observations.  The  agreement  between  the 
two  curves  is  promising  with  differences  being  generally  within  20  W/m2.  There 
are  two  points  to  note  in  this  comparison.  First,  the  maximum  in  the  solar  flux 
occurs  at  approximately  30°  N  for  both  cases.  This  peak  of  net  solar  flux  is  due 
to  the  combined  effects  of  large  insolation  in  the  Northern  Hemisphere  and  small 
cloudiness  in  the  same  regions.  Second,  a  dip  is  found  in  each  of  the  curves  at 
10°  N,  i.e.,  in  the  ITCZ,  in  which  a  maximum  total  cloudiness  is  located.  It  is  also 
noted  that  solar  flux  is  zero  within  the  Antarctic  Circle  because  the  sun  does  not 
rise  above  the  horizon  at  the  region  in  July. 

Comparing  Figures  5.10  and  5.11,  the  model-simulated  OLR  and  net  solar  flux 
at  TOA  are  generally  consistent.  For  instance,  both  the  simulated  OLR  and  net 
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Figure  5.10.  Comparisons  of  the  zonal  mean  OLR  in  the  simulation  and  from  ERB 
data. 
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Figure  5. 11.  Comparisons  of  the  zonal  mean  net  solar  flux  at  TO  A  in  the  simulation 
and  from  EKB  data. 
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solar  flux  exceed  those  from  the  ERB  observations  in  the  tropics. 

5.3. 1.3  Zonal  mean  cloud  LWP 

Data  on  cloud  LWP  are  very  limited.  There  is  no  observational  LWP  data 
on  a  daily  basis  with  which  to  compare,  nor  are  the  zonal  averaged  values  of 
LWP  available  for  July  of  1979.  Since  the  time  domain  in  this  study  is  the  very 
beginning  of  July,  we  have  compared  the  predicted  zonal  mean  LWP  with  the 
monthly  averaged  zonal  mean  data  from  SMMR  for  June,  1979.  Figure  5.12  shows 
the  zonal  mean  LWP  in  the  simulation  and  from  SMMR  measurements  taken  from 
Prabhakara  and  Short  (1984).  Only  the  LWPs  zonally  averaged  over  the  oceans  are 
shown  in  the  figure  for  both  simulation  and  observations.  The  zonal  mean  values 
of  simulated  cloud  LWP  compare  quite  well  with  SMMR  observations.  Differences 
between  the  two  results  are  within  5  mg/cm2  which  is  the  accuracy  of  SMMR  data 
suggested  by  Prabhakara  and  Short  (1984). 

The  maximum  LWP  occurring  in  the  tropics  is  related  to  the  ITCZ.  The  two 
well-defined  local  maxima  of  LWP  at  the  mid-latitudes  in  both  hemispheres  are 
associated  with  the  mid-latitude  storm  tracks.  Minimum  of  LWP  are  found  in 
the  regions  of  subtropical  highs  with  prevailing  downward  motions.  The  SMMR 
data  has  maximum  values  of  about  14,  4.7,  and  4.0  mg/cm2  in  the  tropics  and 
mid-latitudes  of  Southern  and  Northern  Hemispheres,  respectively,  while  the  cor¬ 
responding  values  in  the  simulation  are  13,  5.3,  and  6.5  mg/cm2,  respectively.  The 
ratio  of  maximum  LWP  in  the  tropics  to  those  in  the  middle  latitudes  is  about 
3  to  1.  However,  this  ratio  is  only  about  2  to  1  in  the  model  simulation.  This 
indicates  that  the  ITCZ  peak  in  LWP  is  less  pronounced  in  the  simulation  than  in 
observations.  Obviously,  the  model  underestimates  the  cloud  LWP  in  the  tropics, 
where  a  large  portion  of  LWC  observed  by  the  SMMR  is  condensed  in  convective 
clouds  which  are  not  simulated  in  this  model. 
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Figure  5.12.  Comparisons  of  the  zonal  mean  LWP  in  the  simulation  and  from 
SMMR  data. 


85 


5.3.2  Geographical  distribution  of  cloud  simulations 

The  predicted  cloud  cover,  earth’s  radiation  budget,  LWP,  and  precipitation  in 
the  following  figures  are  presented  in  horizontal  cross-sections  to  demonstrate  the 
geographical  distribution  of  cloud  fields.  Cloud  cover  and  earth’s  radiation  budget 
are  illustrated  by  the  shades  of  gray  with  the  intensity  in  the  picture  proportional 
to  some  ranges  of  data  values. 

5. 3. 2.1  Cloud  cover 

Figures  5.13  and  5.14  show,  respectively,  the  geographical  distribution  of  cloud 
cover  in  the  observation  and  simulation,  with  the  corresponding  gray  values  given 
in  the  lower  right-hand  corner.  The  brighter  the  cloud  picture,  the  larger  the  cloud 
cover  is. 

The  major  distinctive  features  of  the  geographical  distribution  of  cloud  cover 
depicted  in  the  observation  as  well  as  simulation  include: 

1.  A  pronounced  bright  band  of  cloudiness  extending  along  ITCZ  from  150°  W 
to  90°  E,  where  the  cloud  band  splits  into  two  bands.  The  north  cloud  band 
spreads  widely  northeast  into  southern  Asia  and  the  northwestern  Pacific;  the 
south  band  expands  southeast  to  the  southwestern  Pacific.  The  north  band 
is  associated  with  the  seasonal  monsoon  and  the  south  baud  is  related  to  the 
Southern  Pacific  Convergence  Zone  (SPCZ).  Although  the  ITCZ  cloud  band 
is  less  pronounced  in  the  simulation,  the  model  well  predicts  the  bright  cloud 
bands  extending  from  the  central  Pacific  to  the  central  Atlantic,  and  in  central 
Africa. 

2.  Areas  of  small  cloudiness,  i.e.,  a  dark  cloud  picture,  associated  with  subtropical 
highs  occurrng  at  both  north  and  south  of  the  ITCZ  .  For  example,  the 
small  cloud  covers  over  the  southern  Africa,  Europe  and  Middle-East  Asia. 
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Figure  5.14.  The  geographical  distribution  of  total  cloud  cover  in  the  simulation. 
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and  North  and  South  America  are  clearly  identified  in  both  observation  and 
simulation.  The  oceanic  subtropical  highs  in  the  Atlantic  and  Pacific  Oceans 
are  also  evident  by  the  small  cloudiness  over  the  regions. 

3.  The  large  bright  cloud  bands  in  the  mid-latitudes  of  both  hemispheres  as¬ 
sociated  with  the  mid-latitude  storm  tracks.  In  both  the  simulation  and 
observation,  bright  cloud  bands  are  found  over  the  east  coast  of  South  America, 
southern  and  northern  Pacific,  and  the  Greenland  Sea.  The  cloudiness  over  the 
east  coast  of  North  America  is  larger  in  the  simulation  than  in  the  observation. 

4.  The  under-  and  over-estimation  of  cloud  cover  by  3DNEPH  in  the  arctic 
and  antarctic  regions,  respectively,  for  the  summer  season  of  the  Northern 
Hemisphere.  Comparing  Figures  5.13  and  5.14,  it  is  obvious  that  the  cloud 
picture  over  the  antarctic  area  is  brighter  in  the  observation  than  in  the 
simulation,  and  vice  versa  for  the  arctic  region. 

5. 3. 2. 2  Earth’s  radiation  budget 

The  earth’s  radiative  flux  at  TOA  reflects  many  of  the  features  in  the  simulated 
cloudiness  discussed  in  the  previous  section.  It  is  the  radiative  effects  of  clouds 
rather  than  the  cloud  cover  alone  that  ultimately  determines  the  model’s  simulation 
of  the  radiation  budget.  Therefore,  comparisons  of  the  simulated  earth’s  radiation 
budget  with  ERB  data  may  provide  valuable  guidance  on  the  cloud  radiative 
properties  of  the  LWC/IWC  simulated  in  the  model. 

Figures  5.15  and  5.16  show  the  outgoing  longwave  radiation  in  the  observation 
and  simulation,  respectively.  Figures  5.17  and  5.18  are  the  same  as  Figures  5.15 
and  5.16,  respectively,  except  for  net  solar  downward  flux  at  TOA.  The  intensity  of 
the  gray  shades,  given  in  the  lower  right-hand  corner  of  each  figure,  is  defined  such 
that  the  brighter  the  picture,  the  smaller  the  radiative  flux. 
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Figure  5.15.  The  geographical  distribution  of  OLR  from  ERB  data. 
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Figure  5.16.  The  geographical  distribution  of  OLR  in  the  model  simulation. 
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Figure  5.18.  The  geographical  distribution  of  simulated  net  solar  flux  at  TO  A. 
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In  all  four  figures,  small  flux  values  from  the  90°  S  to  about  40°  S  are  due  to 
relatively  low  solar  radiation  in  the  winter  hemisphere.  Therefore,  very  few  flux 
features  over  these  regions  are  evident  in  Figures  5.15  -  5.18.  However,  north  of  the 
40°  S,  more  flux  features  are  evident  in  all  of  the  figures  because  of  relatively  large 
solar  zenith  angle  and  more  variability  in  cloud  cover. 

The  earth’s  radiation  budget  in  the  simulation  shows  good  agreement  with  that 
in  the  observation.  There  are  two  distinctive  features  on  these  maps.  First,  the 
well-defined  ITCZ  and  monsoon  circulation  are  characterized  by  a  bright  picture, 
i.e.,  small  values  of  OLR  and  net  solar  flux.  They  are  located  over  the  equatorial 
western  and  eastern  Pacific,  the  Atlantic,  and  central  Africa,  as  well  as  the  monsoon 
area  over  the  India  and  Asia.  As  mentioned  previously,  the  cloudiness  in  these 
regions  is  underestimated  in  the  simulation.  As  would  be  expected,  the  flux  values 
in  these  regions  are  overestimated  by  the  model.  Second,  the  almost  continuous 
belts  of  dark  picture  are  associated  with  the  subtropical  highs  located  to  the  south 
and  north  of  the  ITCZ.  In  these  regions,  large  flux  values  for  both  OLR  and  net  solar 
flux  are  due  to  the  small  cloudiness  dominating  the  areas.  Small  flux  values  arc  also 
evident  in  the  areas  of  the  mid-latitude  storm  tracks  in  the  Northern  Hemisphere. 
However,  the  simulated  values  of  OLR  and  net  solar  flux  in  the  Greenland  Sea, 
near  70°  N  and  0°  E,  appear  to  be  too  small  compared  with  the  observation. 

5. 3. 2. 3  Liquid  water  path 

Figures  5.19  and  5.20  illustrate  the  geographical  distribution  of  LWP  in  the 
observation  and  simulation,  respectively.  The  observation,  taken  from  Prabhakara 
and  Short(1984),  is  a  monthly  averaged  map  for  July,  1979.  As  was  mentioned 
previously,  the  SMMR  observations  were  limited  to  the  oceans  equatorward  of  60° 
latitude  because  of  a  nearly  constant  emissivity  over  the  ice-free  sea  surface.  As 
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would  be  expected,  Figure  5.19  shows  no  SMMR  results  on  LWP  over  land  and  at 
high  latitudes. 

In  both  of  the  Figures,  large  amounts  of  liquid  water  placed  in  the  ITCZ  and 
Asian  monsoon  regions  are  clearly  evident,  while  in  the  subtropical  highs  at  both  of 
hemispheres  small  amounts  of  liquid  water  are  found.  In  the  central  and  southern 
Pacific,  the  model  simulates  large  values  of  LWP,  which  are  also  evident  in  the 
observation.  In  both  the  simulation  and  observation,  an  abundance  of  liquid  water 
is  found  over  the  Indian  Ocean,  and  southeast  and  northeast  Asia.  In  addition, 
the  well-defined  observed  local  maxima  of  LWP  near  the  east  coasts  of  North, 
Central,  and  South  America  are  clearly  simulated  by  the  model.  Localized  LWP 
minima  observed  west  of  Africa,  and  South  and  North  America  are  also  shown  in  '  he 
simulation.  These  minima  are  associated  with  large-scale  atmospheric  subsidences; 
a  mechanism  which  acts  to  prevent  water  vapor  from  condensing,  greatly  reducing 
cloud  LWC  in  these  regions. 

The  large  amounts  of  liquid  water  associated  with  the  mid-latitude  storm  tracks 
near  50°  S(N)  -  60°  S(N)  are  not  shown  in  the  SMMR  observation  because  of  high 
latitudes  where  sea  surfaces  are  covered  with  ice.  In  addition,  the  model  simulates 
large  amounts  of  liquid  water  over  land  at  Central  America,  central  Africa,  India  and 
southeast  China.  These  results  are  not  reported  by  the  SMMR  observations  due  to 
the  inhomogeneous  emissivity  over  the  land  surface.  From  the  above  comparisons, 
it  appears  that  the  geographic  distribution  of  LWP  is  well  simulated. 

5.3. 2. 4  Precipitation 

Precipitation  is  initiated  by  cloud  microphysical  processes  such  as  collision  and 
coalescence,  through  which  small  cloud  droplets  merge  into  large  raindrops.  How¬ 
ever,  the  model’s  large-scale  precipitation  distribution,  when  averaged  over  ex¬ 
tended  periods,  is  determined  by  the  circulation  in  the  model,  rather  than  by 
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Figure  5.20.  The  geographical  distribution  of  LWP  in  the  simulation. 


the  cloud  microphysical  processes  themselves  (Smith’s,  1990).  Moreover,  observed 
precipitation  is  largely  generated  by  convective  clouds,  which  are  not  included  in 
the  model  simulation.  Due  to  the  simplicity  of  the  prescribed  wind  fields  and 
the  exclusion  of  convective  clouds  in  this  study,  it  is  not  our  intent  to  verify 
quantitatively  the  simulated  precipitation  at  this  stage.  Instead,  the  geographical 
distribution  of  simulated  precipitation  is  shown  in  this  section. 

Figure  5.21  shows  the  geographical  distribution  of  simulated  precipitation.  In 
the  figure,  a  noticeable  rain  belt  is  located  along  the  ITCZ.  The  model  well  simulates 
the  precipitation  maxima  of  both  the  Atlantic  and  the  Pacific  branches  of  the  ITCZ, 
as  well  as  that  in  the  South  Pacific  Convergence  Zone.  The  precipitation  minima 
are  simulated  in  the  regions  of  subtropical  highs,  located  south  and  north  of  the 
ITCZ.  Little  precipitation  is  simulated  over  the  deserts,  such  as  the  Sahara  and  the 
Gobi.  The  wide-spread  rain  bands  associated  with  the  mid-latitude  storm  tracks 
are  shown  in  the  simulation.  Large  amounts  of  simulated  precipitation  is  shown  on 
the  east  coast  of  America,  in  contrast  to  small  amounts  of  precipitation  on  the  west 
coast.  A  noticeable  amount  of  precipitation  is  simulated  over  central  Africa,  while 
little  precipitate  can  be  found  over  west  Africa. 

5.3.3  Cloud  ice  content 

Since  observational  data  on  cloud  IWC  is  not  available  on  the  global  basis,  we 
qualitatively  compare  the  model-predicted  cloud  IWC  with  aircraft  measurements 
in  cirrus,  as  well  as  with  other  model  results.  Based  on  the  aircraft  measurements  in 
mid-latitude  cirrus  over  the  United  States,  Heymsfield  and  Platt  (1984)  derived  the 
dependence  of  cirrus  cloud  properties,  including  cloud  IWC,  on  the  temperature. 
According  to  their  report,  averaged  cloud  IWC  varied  from  a  few  milligrams  per 
cubic  meter  at  temperatures  below  —40°  C  to  about  20  mg/m3  above  —40°  C. 
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Similar  magnitudes  of  cloud  IWC  were  also  simulated  in  the  cirrus  cloud  model  by 
Starr  and  Cox  (1985). 

The  zonal  average  of  simulated  cloud  IWC,  averaged  over  July  3-4,  1979,  is 
shown  in  Figure  5.22.  The  magnitude  of  predicted  cloud  IWC  is  on  the  order  of  10 
n  •  'm3,  which  is  consistent  with  the  aforementioned  model  simulation  and  aircraft 
obsei . .  ‘  ,,al  values  in  cirrus. 

Comparisons  between  Figures  5.22  and  4.1  indicate  the  maxima  of  cloud  IWC 
locate  along  —  15°  C,  in  which  the  maximum  Bergeron-Findeison’s  processes  occur. 
The  simulated  cloud  IWC  decreases  from  about  10  mg/m3  near  —  15°  C  to  a 
few  milligrams  per  cubic  meter  at  temperatures  below  —  40°  C.  The  variations 
of  cloud  IWC  with  respect  to  the  temperature  are  similar  to  those  suggested  by 
Heymsfield  and  Platt  (1984).  It  is  noted  that  the  simulated  cloud  ice  disappears 
at  temperatures  warmer  than  the  freezing  point,  since  the  melting  process  converts 
ice  crystals  into  cloud  liquid  water.  Ice  crystals  and  liquid  water  may  co-exist  at 
the  region  above  the  freezing  level. 

The  latitudinal  distribution  of  cloud  IWC  appears  realistic.  In  the  antarctic  area, 
Southern  Hemisphere  wintertime,  all  low,  middle,  and  high  clouds  are  ice  clouds, 
while  in  the  arctic  region,  Northern  Hemisphere  summertime,  only  the  middle  and 
high  clouds  are  ice  clouds.  In  the  tropics,  only  high  clouds  are  dominated  by  ice 
crystals. 

Based  on  the  above  discussions,  the  simulated  cloud  IWC  appears  to  depend 
more  on  temperature  than  on  vertical  velocity,  which  is  important  to  the  conden¬ 
sational  process.  The  predicted  cloud  IWC  appears  to  agree  very  well  with  that  of 
the  aircraft  observations  in  cirrus  clouds. 

5.4  Discussion 

The  comparisons  shown  in  the  previous  sections  indicate  the  developed  large- 

scale  cloud  model  is  able  to  simulate  rather  realistically  zonal  means  and  geo- 

100 


graphical  distributions  of  cloud  fields,  including  cloud  cover,  cloud  LWC/ 1 WC  and 
precipitation.  By  interactively  providing  the  cloud  cover,  cloud  LWC  and  IWC 
to  the  radiative  transfer  calculations,  the  earth’s  radiation  budget,  namely  OLR 
and  net  solar  flux  at  TOA,  is  also  well  simulated.  More  importantly,  the  previous 
comparisons  demonstrate  the  consistency  between  the  simulated  OLR  and  net  solar 
flux  at  the  TOA.  For  instances,  both  of  them  are  small  over  the  regions  of  ITCZ 
and  monsoon,  and  large  at  the  subtropical  highs. 

While  the  comparisons  between  the  simulations  and  observations  attempt  to 
verify  the  model  prediction,  intercomparisons  among  the  simulated  fields,  including 
cloud  cover  (Fig.  5.14),  OLR  (Fig.  5.16),  net  solar  flux  (Fig.  5.18),  LWP  (Fig.  5.20) 
and  precipitation  (Fig.  5.21)  may  render  us  some  physical  insight  into  the  relation¬ 
ship  between  the  large-scale  circulations  and  cloud  formation. 

The  model  simulation  shows  that  clouds  are  strongly  tied  to  the  large-scale 
circulations.  For  instances,  the  influence  of  the  Asian  monsoon  circulation  and 
ITCZ  on  cloud  formation  are  prominent.  These  regions  are  associated  with  con¬ 
vergences  of  low-level  winds  (i.e.,  upward  motions),  and  deep  moist  tropospheres. 
The  low-level  wind  convergences  result  in  upward  motions  and  cooling  expansions 
of  the  air  parcels.  At  the  same  time,  the  deep  moist  tropospheres  provide  the  water 
vapor  that  allows  condensation  to  occur  in  the  air  parcels.  Consequently,  large 
amounts  of  cloud  form  in  these  regions,  as  evident  in  Figures  5.14,  5.20  and  5.21. 
Observing  the  model-simulated  results  in  these  figures  for  the  ITCZ  and  monsoon 
areas,  we  find  model  consistently  simulates  large  amounts  of  cloud  cover,  LWP, 
and  precipitation  in  these  regions.  Large  amounts  of  cloud  cover  coincide  with 
large  amounts  of  cloud  LWC,  which  generate  large  precipitation  rates  by  virtue  of 
autoconversion.  Similarly,  large  amounts  of  cloud  are  simulated  in  the  regions 
of  the  mid-latitude  storm  tracks,  which  are  characterized  by  strong  baroclinic 
instability  in  the  atmosphere.  These  baroclinic  disturbances  are  accompanied  by 
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cloudiness  maxima,  abundance  of  liquid  water,  and  noticeable  rain  bands  in  the 
middle  latitudes  of  both  hemisphere. 

At  the  subtropical  highs  of  both  hemispheres,  small  amounts  of  cloud  cover, 
LWP  and  precipitation  in  the  simulation  are  associated  with  large-scale  atmospheric 
subsidence,  which  is  related  to  the  branches  of  downward  motion  in  the  Hadley  cir¬ 
culations.  The  downward  motion  tend  to  increase  the  atmospheric  temperature  and 
result  in  the  decrease  of  relative  humidity,  which  acts  to  prevent  water  vapor  from 
condensing.  As  a  result,  small  cloudiness,  LWP  minima,  and  little  precipitation 
are  simulated  in  these  regions.  It  is  also  evident  from  the  simulation  that  the  east 
coast  of  America  is  favorable  for  cloud  formation,  while  the  west  coast  is  not.  Large 
amounts  of  cloud  cover,  liquid  water,  and  precipitation  are  simulated  on  the  east 
coast  of  America,  but  little  cloud  forms  in  the  west  coast. 

One  of  the  main  problems  in  the  verification  of  cloud  simulation  against  satellite 
data  is  associated  with  the  comparisons  of  vertically  integrated  liquid  water  (LWP). 
The  LWP  is  given  by  the  cloud  LWC  and  model  thickness  as  follow 

LWP  =  LWC  x  Az 

where  LWC  is  the  cloud  LWC  denoted  earlier  and  Az  is  the  model  vertical  grid 
thickness.  In  the  context  of  GCMs,  clouds  are  sub-grid  scale  in  the  horizontal,  and 
more  seriously,  in  the  vertical.  However,  in  all  of  the  GCM  simulations,  clouds  are 
assumed  to  fill  the  entire  vertical  grid  layer  whenever  clouds  form.  Consequently, 
the  LWPs  in  the  simulation  are  overestimated  due  to  the  exaggeration  of  cloud 
vertical  extent  in  the  model. 

For  example,  low  cloud  base  and  middle  cloud  top  in  the  model  are  defined 
at  925  mb  and  450  mb,  respectively.  Thus  the  thickness  of  the  low  and  middle 
clouds,  in  which  most  cloud  LWCs  are  located,  may  be  roughly  estimated  as  5  km. 
The  magnitude  of  cloud  LWC,  based  on  the  aircraft  observations  in  stratus  clouds 
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(Kriebel,  1989),  may  be  estimated  as  0.1  g/m3,  i.e.,  10“4  mg/cm3,  which  is  on  the 
order  of  the  stratiform  cloud  LWC  adopted  in  Liou  and  Wittman  (1977).  The  LWP, 
given  by  the  product  of  10-4  mg/cm3  and  5  x  105  cm,  is  as  large  as  50mg/cm2, 
which  is  about  ten  times  larger  than  that  derived  from  the  SMMR  observations  in 
middle  latitudes.  Thus,  in  order  to  obtain  reasonable  LWP  with  the  coarse  model 
thickness,  the  cloud  LWC  has  to  be  smaller  than  it  is  actually  observed  to  be  in 
the  aircraft  observations. 

With  the  unavoidable  overestimation  of  cloud  thickness,  the  decision  needs  to 
be  made  whether  the  model  attempts  to  simulate  reasonable  values  in  LWP  or  in 
cloud  LWC.  The  motive  of  forecasting  cloud  LWC  is  to  calculate  interactively  cloud 
radiative  properties.  Furthermore,  since  it  is  the  effects  of  vertically-integrated 
cloud  liquid  water  content  (i.e.,  LWP),  rather  than  cloud  LWC  alone  that  ultimately 
determines  the  cloud  radiative  properties,  we  have  decided  that  the  reasonable 
amount  of  LWP  is  more  important  than  that  of  cloud  LWC.  In  addition,  there  is 
no  observational  cloud  LWC  that  can  be  utilized  as  verification  data  on  the  global 
basis.  The  model-simulated  cloud  LWCs  are  verified  against  SMMR  observations 
in  terms  of  LWP. 

Having  recognized  the  importance  of  LWP,  we  have  adopted  a  large  value  for  the 
autoconversion  rate  to  reduce  the  amount  of  cloud  LWC  in  the  model,  so  that  the 
model-simulated  LWP  may  be  comparable  with  the  SMMR  observations.  Unless 
specified,  the  autoconversion  rate  is  adopted  as  Co  =  5  x  10-4s“1,  which  is  5  times 
larger  than  that  used  in  Sundqvist  (1978). 
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Section  6 


SENSITIVITY  STUDIES 


In  order  to  examine  the  interactions  and  feedbacks  of  radiative  heating  fields  and 
ice  phase  cloud  microphysics  on  the  large-scale  cloud  simulations,  a  set  of  numerical 
experiments  has  been  carried  out  with  and  without  the  inclusion  of  radiation  and 
ice  phase  processes.  A  brief  description  of  the  experimental  runs  follows: 

•  CTRL:  Control  run  which  includes  both  radiation  and  ice-phase  processes. 

•  EXP1:  As  in  CTRL  but  without  the  radiation. 

•  EXP2:  As  in  EXPl  but  without  the  ice-phase  processes. 

All  of  these  control  and  experimental  runs  have  the  same  initial  conditions,  and 
wind  fields,  and  their  results  are  averaged  over  the  same  periods,  allowing  inter¬ 
comparisons  to  be  made  between  the  control  run  and  the  experimental  runs. 

6.1  Control  Run 

The  control  run  is  the  same  run  as  that  described  in  the  previous  section.  Its 
results  were  shown  in  the  last  section  on  the  horizontal  cross  sections.  In  this 
section,  these  results  are  demonstrated  in  meridional  cross  section  to  show  the 
vertical  structure  of  cloud  fields. 

Figures  6.1  -  6.4  show  the  zonal  average  in  temperature,  cloud  cover,  cloud  LWC 
and  cloud  IWC  predicted  and  averaged  over  July  3-4,  1979,  in  the  control  run. 
The  profiles  of  predicted  cloud  cover,  shown  in  Figure  6.2,  are  related  to  those  of 
vertical  velocity  depicted  in  Figure  4.5. 
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he  zonal  mean  field  of  cloud  cover  in  the  control  run. 
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f  cloud  LWC  in  the  control 


Three  large  amounts  of  cloud  cover,  located  in  the  tropics  and  at  the  mid-high 
latitudes  at  both  of  hemispheres,  are  due  to  the  large  upward  motions  in  these 
regions,  while  two  relatively  clear-sky  areas  at  the  subtropical  highs  in  both  hemi¬ 
spheres  are  related  to  the  downward  motions  in  the  regions.  In  addition,  cloud 
cover  in  the  southern  hemispheric  subtropics  is  smaller  than  that  in  the  North¬ 
ern  Hemisphere  because  of  the  branch  of  strong  sinking  motion  of  the  Southern 
Hemisphere  Hadley  cell. 

The  pattern  of  simulated  cloud  LWC  in  Figure  6.3  is  consistent  with  that  of 
cloud  cover  in  Figure  6.2;  i.e.,  cloud  LWC  exists  wherever  there  is  cloud.  The 
production  of  cloud  LWC  in  the  simulation  depends  on  the  vertical  velocity  and 
amount  of  moisture  in  the  region.  If  the  moisture  supply  is  the  same,  the  stronger 
the  upward  motion,  the  more  cloud  LWC  is  generated  due  to  a  large  condensational 
rate  resulting  from  a  large  vertical  velocity.  In  the  vertical,  large  amounts  of  cloud 
LWC  are  found  in  low  and  middle  clouds  because  of  abundant  moisture  available  to 
condensation  in  low  and  middle  atmosphere.  Similarly,  the  maximum  cloud  LWC 
is  located  at  the  tropics  and,  two  local  maxima  exist  at  the  middle  latitudes  in  both 
hemispheres. 

The  three  large  amounts  of  cloud  LWC,  located  at  the  tropics  and  mid-latitude 
storm  tracks  in  both  hemispheres,  correspond  to  the  large  cloudiness  in  these 
regions.  The  small  amounts  of  cloud  LWC  associated  with  the  subtropical  highs 
are  due  to  the  prevailing  downward  motions  in  the  regions.  In  the  antarctic  area, 
clouds  contain  small  amounts  of  LWC  because  of  the  low  temperatures  in  the  area. 

The  zonal  mean  of  cloud  IWC,  simulated  in  the  control  run,  is  shown  in  Fig¬ 
ure  6.4,  which  is  identical  to  Figure  5.22.  Comparing  of  Figures  6.4  and  6.3,  in 
the  regions  of  high  latitude  or  high  altitude,  where  temperatures  are  low,  clouds 
are  dominated  by  the  ice  crystals,  while  in  the  tropics  and  low  latitudes,  where 
temperatures  are  higher,  cloud  LWC  prevails  in  the  clouds. 
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6.2  Experiment  On  Inclusion  Of  Radiation 

Figures  6.5  -  6.12  show  the  results  from  the  experiment  with  and  without  the 
inclusion  of  radiation.  These  figures  indicate  the  differences  obtained  by  subtract¬ 
ing  the  results  in  EXP1,  i.e.,  without  radiation,  from  those  in  CTRL,  i.e.,  with 
radiation. 

The  differences  in  the  solar  radiative  heating  rates  are  illustrated  in  Figure  6.5, 
which  show  all  the  positive  values.  The  results  indicate  that  solar  radiation  basically 
heats  the  atmosphere.  As  would  be  expected,  the  solar  heating  rates  in  July 
in  the  Northern  Hemisphere  are  much  larger  than  in  the  Southern  Hemisphere. 
Furthermore,  there  is  no  solar  radiation  within  the  Antarctic  Circle  because  the 
sun  does  not  rise  above  the  horizon  in  the  region  in  July. 

There  are  three  heating  rate  maxima  with  values  of  about  2°  K/day  located  at 
the  tops  of  middle  cloud  in  the  Northern  Hemisphere.  The  large  values  of  solar 
heating  rates  within  the  Arctic  Circle  are  due  to  a  longer  solar  day  in  the  region 
in  July,  while  large  solar  heating  rates  in  the  tropics  are  due  to  large  solar  zenith 
angles  during  the  daytime.  The  large  solar  heating  located  in  the  upper  atmosphere 
is  associated  with  the  strong  O3  absorption  at  high  altitudes. 

Figure  6.6  is  the  same  as  Figure  6.5,  except  for  the  IR  radiative  cooling  rates. 
The  most  distinctive  feature  in  the  figure  is  the  strong  IR  cooling  rates  occurring 
at  the  middle  cloud  tops.  This  is  because  high  cloud  is  considered  nonblack  in 
the  radiation  calculations,  while  middle  cloud  is  assumed  to  be  a  blackbody.  As 
a  result,  the  upward  IR  flux  emitted  by  middle  cloud  is  much  stronger  than  the 
downward  flux  emitted  by  high  cloud.  Consequently,  strong  IR  radiative  cooling  is 
generated  at  the  top  of  middle  cloud. 

As  evident  from  Figure  6.6,  the  inclusion  of  IR  radiation  basically  cools  the 
atmosphere.  However,  if  low  cloudiness  is  large  enough  to  emit  a  strong  downward 
IR  flux  out  of  the  cloud  base,  weak  IR  heating  may  be  produced  beneath  the  cloud 
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IR  HEATING  RATE  DIFFERENCE  (degree  C/day) 


Figure  6.6.  Zonally  averaged  differences  in  IR  radiative  heating  rates.  The 
differences  are  obtained  by  subtracting  the  results  in  EXP1  from  those  in  CTRL. 
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base.  This  explains  the  two  regions  of  weak  IR  heating  at  the  low  cloud  bottoms 
at  60°  S  and  in  the  tropics.  Two  other  IR  heating  rates,  located  in  the  upper 
atmosphere  in  the  tropics  and  in  the  lower  atmosphere  in  the  Arctic,  are  associated 
with  the  areas  of  high  ozone  concentration  and  temperature  inversion  in  the  regions, 
respectively. 

Figure  6.7  is  the  same  as  Figure  6.5  except  for  the  net  radiative  heating  rates 
which  are  the  combined  effects  of  solar  and  IR  radiation.  The  most  important 
conclusion  drawn  from  Figure  6.7  is  that  the  inclusion  of  the  radiative  processes 
basically  decreases  the  atmospheric  temperature.  As  would  be  expected,  the  net 
radiation  shows  relatively  large  cooling  at  middle  cloud  tops  and  small  cooling  at 
cloud  bottoms.  However,  there  are  weak  heating  rates  in  the  low  clouds  in  the 
tropics  and  Southern  Hemisphere  due  to  large  low  cloud  cover  over  those  regions. 
The  net  radiative  cooling  is  larger  in  high  clouds  than  in  low  clouds. 

The  strongest  radiative  cooling  at  the  middle  cloud  top,  located  within  the 
Antarctic  Circle,  is  due  to  the  absence  of  solar  radiation  in  the  winter  season. 
Relatively  large  cooling  at  the  low  clouds  within  the  arctic  area  is  associated  with 
the  temperature  inversion  in  the  region,  while  the  large  cooling  near  the  surface 
layers  in  the  equatorial  region  is  due  to  the  contribution  of  water  vapor  conlinnum 
absorption. 

It  is  noted  that  noticeable  net  radiative  heating  is  located  in  the  upper  atmo¬ 
sphere  in  the  Northern  Hemisphere  because  of  the  ozone  solar  absorption  in  the 
region.  The  weak  heating  in  middle  clouds  within  the  Arctic  Circle  is  due  to 
absorption  of  solar  radiation  during  the  long  arctic  solar  day. 

Figure  6.8  is  the  same  as  Figure  6.5  except  for  the  temperatuie.  In  Figure  6.8,  the 
major  features  are  the  strong  temperature  decreases  at  the  middle  cloud  tops  and 
the  relatively  small  temperature  reductions  in  the  low  clouds.  This  is  because  the 
radiative  cooling  is  much  stronger  at  a  middle  cloud  top  than  a  low  cloud  bottom. 
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Figure  6.7.  Zonally  averaged  differences  in  net  radiative  heating  rates.  The 
differences  are  obtained  by  subtracting  the  results  in  EXP1  from  those  in  CTRL. 
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Figure  6.8.  Zonally  averaged  differences  in  temperature.  The  differences  are 
obtained  by  subtracting  the  results  in  EXP1  from  those  in  CTRL. 
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The  strongest  temperature  reduction  occurs  at  the  middle  cloud  top  near  60°  S,  at 
which  the  sun  does  not  rise  above  the  horizon  in  the  winter  season.  Temperature 
decreases  are  relatively  small  in  the  subtropical  regions  where  small  cloud  amounts 
are  located. 

The  inclusion  of  radiation  decreases  the  atmospheric  temperature;  i.e.,  there  is 
a  cooling  effect  on  the  atmosphere.  The  temperature  decrease  is  larger  for  cloudy 
sky  than  clear  sky.  These  changes  in  temperature  are  contributed  by  the  feedbacks 
of  the  net  radiative  cooling  in  Figure  6.7  through  the  thermodynamic  equation. 
These  temperature  changes,  in  turn,  initiate  all  of  the  changes  in  the  cloud  fields 
such  as  cloud  cover,  cloud  LWC,  and  cloud  IWC. 

Figure  6.9  is  the  same  as  Figure  6.5  except  for  the  differences  in  cloud  cover. 
Comparisons  between  Figures  6.9  and  6.5  indicate  that  the  patterns  of  differences 
in  cloud  cover  and  temperature  are  similar.  They  are  related  such  that  decreases  in 
temperature  result  in  increased  cloud  cover,  and  vice  versa.  A  decrease  in  temper¬ 
ature  initiates  two  competitive  effects  on  the  relative  humidity,  which  ultimately 
determines  the  cloud  fraction.  One  effect  is  the  increase  of  condensation,  which 
depletes  more  water  vapor  and  results  in  reducing  relative  humidity.  The  other  is 
the  strong  decrease  of  saturation  vapor  pressure,  which  greatly  increases  relative 
humidity.  It  is  evident  from  the  Figure  6.9  that  clouds  are  dominated  by  the  latter 
effect  which  results  in  the  increase  of  cloud  cover  with  decreasing  temperature. 

The  inclusion  of  radiation  increases  cloud  cover.  The  large  increases  in  cloud 
cover  occur  at  the  middle  cloud  tops,  in  which  large  temperature  decreases  are 
located.  The  increases  in  low  cloudiness  is  relatively  small  due  to  small  radiative 
cooling  associated  with  low  clouds.  However,  there  are  two  large  increases  of  low 
cloud  cover.  One  occurs  in  the  tropics  where  a  large  amount  of  water  vapor  exists, 
while  the  other  is  in  the  arctic  area  where  a  strong  inversion  is  located. 

Figures  6.10,  6.11  and  6.12  are  the  same  as  Figure  6.5  except  for  the  differences 
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Figure  6.9.  Zonally  averaged  differences  in  cioud  cover.  The  differences  are  obtained 
by  subtracting  the  results  in  EXP1  from  those  in  CTRL. 
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averaged  differences  in  cloud  LWC.  The  differences  are 
the  results  in  EXP1  from  those  in  CTRL. 
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Figure  6.11.  Zonally  averaged  differences  in  cloud  IWC.  The  differences  are 
obtained  by  subtracting  the  results  in  EXPl  from  those  in  CTRL. 
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Figure  6.12.  Zonally  averaged  differences  in  total  cloud  water.  The  differences  are 
obtained  by  subtracting  the  results  in  EXP1  from  those  in  CTRL. 
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in  cloud  LWC,  IWC,  and  total  cloud  water  content,  respectively.  As  was  mentioned, 
total  cloud  water  is  the  combination  of  cloud  LWC  and  IWC.  Therefore,  Figure  6.12 
is  the  combined  picture  of  Figures  6.10  and  6.11.  Cloud  water  increases  wherever 
cloud  cover  increases.  As  was  mentioned  previously,  the  decrease  of  cloud  tem¬ 
perature  by  radiative  cooling  increases  the  condensational  rate  and  results  in  the 
increase  of  total  cloud  water.  The  large  increase  of  cloud  water  occurs  at  the  tops 
of  middle  and  low  cloud  in  the  tropics,  where  strong  radiative  cooling  is  located. 

6.3  Experiment  On  Ice  Phase  Processes 

Figures  6.13  -  6.17  show  the  results  from  the  experiment  with  and  without  the 
inclusion  of  ice-phase  processes.  These  figures  depict  the  differences  obtained  by 
subtracting  the  results  in  the  EXP2,  i.e.,  without  ice-phase  processes,  from  those 
in  the  EXP1,  i.e.,  with  ice-phase  processes.  It  is  noted  that  radiative  processes  are 
not  included  \u  these  experimental  runs. 

The  differences  in  temperature  and  cloud  cover  due  to  the  inclusion  of  ice-phase 
processes,  shown  in  Figures  6.13  and  6.14,  are  insignificant.  Figure  6.15  illustrates 
the  differences  in  the  cloud  LWC  between  the  cases  with  and  without  the  inclusion  of 
ice-phase  processes.  The  figure  shows  little  change  in  cloud  LWC  in  the  tropical  low 
and  middle  clouds,  in  which  ice  clouds  are  less  dominant.  The  marked  reductions 
of  cloud  LWC  occur  in  low,  middle,  high  clouds  within  the  Antarctic  Circle,  high 
clouds  in  the  tropics,  and  middle,  and  high  clouds  within  the  Arctic  Circle.  In  these 
regions  clouds  are  dominated  by  ice  crystals,  which  deplete  the  cloud  LWC  by  the 
virtue  of  Bergeron-Findeison’s  processes.  Figure  6.16  is  the  same  as  Figure  6.15 
except  for  the  differences  in  cloud  IWC.  Since  there  is  no  cloud  IWC  when  ice 
phase  processes  are  excluded,  the  differences  in  cloud  IWC  shown  in  Figure  6.16 
are  identical  to  the  cloud  IWCs  shown  in  the  control  run  (see,  Figure  6.4). 

Figure  6.17  is  the  same  as  Figure  6.15  except  for  the  differences  in  total  cloud 

water  content,  i.e.,  cloud  LWC  plus  IWC.  The  differences  in  the  total  cloud  water 

122 


-60.  -30.  0.  30.  60.  90. 

LATITUDE  (DEGREE) 


Zonally  averaged  differences  in  temperature.  The  differences  are 
btracting  the  results  in  EXP2  from  those  in  EXP  1. 
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Figure  6.15.  Zonally  averaged  differences  in  cloud  LWC.  The  differences  are 
obtained  by  subtracting  the  results  in  EXP2  from  those  in  EXP1. 
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content  is  expected  to  be  smaller  than  in  the  cloud  LWC,  since  the  existence  of 
cloud  IWC  compensates  for  the  decreased  cloud  LWC.  However,  Figure  6.17  shows 
that  the  amounts  of  cloud  IWC  not  only  partially  offset  the  loss  of  cloud  LWC,  but 
overcompensate  the  depleted  cloud  LWC  in  the  regions  where  ice  cloud  prevails. 

The  maximum  increases  in  total  cloud  water  content  occur  near  60°  S,  40°  N, 

N 

and  60°  —  85°  N.  The  increases  of  total  cloud  water  when  the  ice  phase  is  included 
indicate  the  existing  ice  crystals  in  these  regions  are  not  generated  by  the  Bergeron- 
Findeison’s  processes  alone.  In  fact,  large  amounts  of  the  increased  cloud  IW C  are 
due  to  advection.  For  examples,  comparisons  of  Figures  6.17  and  6.4  indicate 
the  maximum  increase  of  total  cloud  water  at  40°  N  corresponds  to  the  minimum 
cloud  IWC  in  the  region,  where  the  Bergeron- Findeison’s  processes  are  relatively 
weak.  In  addition,  comparisons  between  Figures  6.17  and  4.5  suggest  the  maximum 
increases  of  total  cloud  water  are  associated  with  upward  motions,  i.e.,  convergence, 
near  60°  S,  40°  N,  65°  N  and  80°  N. 

From  the  above  discussions,  it  is  evident  that  the  reduction  of  cloud  LWC 
by  the  Bergeron-Findeison’s  processes  is  less  effective  than  the  increase  of  cloud 
IWC  due  to  advection.  This  is  because  the  simulated  cloud  LWC  is  too  small  to 
provide  effective  Bergeron-Findeison’s  processes  to  overcome  the  advection  effect. 
For  the  reasons  described  in  Section  5.4,  the  autoconversion  rate  adopted  in  EXP1 
and  EXP2  are  five  times  larger  than  that  used  in  Sundqvist  (1978).  The  large 
autoconversion  rate  (Co  =  5  x  10-4)  adopted  in  this  study  is  responsible  for  the 
small  cloud  LWC  in  the  simulation.  As  was  mentioned  previously,  the  simulated 
cloud  LWC  is  subject  to  the  adopted  value  of  the  autoconversion  rate;  the  smaller 
the  autoconversion  rate,  the  larger  the  simulated  cloud  LWC. 

In  order  to  further  examine  the  effects  of  Bergeron-Findeison’s  processes  on 
cloud  formation,  we  have  undertaken  two  more  experimental  runs  with  the  same 
autoconversion  rate  as  that  used  in  Sundqvist  (1978),  i.e.,  C0  =  10"4.  Thus 
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the  simulated  cloud  LWC  is  large  enough  to  provide  strong  Bergeron-Findeison’s 
processes.  In  the  following  two  experimental  runs,  the  first  run  is  the  same  as  EXPl 
with  the  smaller  value  of  the  autoconversion  rate,  i.e.,  Co  —  10-4;  the  second  run 
is  the  same  as  the  first  run  but  without  the  ice  phase  processes. 

Figures  6.18  -  6.22  are  the  same  as  Figures  6.13  -  6.17,  except  for  the  smaller 
value  of  the  autoconversion  rate.  The  temperature  changes  in  Figure  6.18  are 
similar  to  those  in  Figure  6.13.  However  the  temperature  changes  in  the  former  are 
larger  than  in  the  latter.  This  indicates  that  the  Bergeron-Findeison’s  processes 
are  stronger  with  a  smaller  autoconversion  rate  than  they  are  with  larger  rate.  The 
differences  in  cloud  cover,  shown  in  Figure  6.19,  are  insignificant.  Figure  6.20  shows 
a  very  similar  pattern  as  Figure  6.15.  Due  to  a  smaller  value  of  the  autoconversion 
rate  in  Figure  6.20,  the  values  of  the  differences  for  cloud  LWC  are  about  5  times 
larger  in  Figure  6.20  than  they  are  in  Figure  6.15.  The  large  differences  in  cloud 
LWC  in  Figure  6.20  indicate  strong  Bergeron-Findeison’s  processes  happen  when 
there  is  more  cloud  LWC.  However,  the  values  of  cloud  IWC  in  Figure  6.21  are  close 
to  those  in  Figure  6.16,  even  though  the  Bergeron-Findeison’s  processes  are  much 
stronger  in  the  former  than  in  the  latter.  This  indicates  the  downward  ice  flux  is 
larger  in  the  former  than  in  the  latter. 

Figure  6.22  is  the  same  as  Figure  6.20  except  for  the  differences  in  total  cloud 
water  content.  Figure  6.22  shows  a  very  similar  pattern  as  Figure  6.20.  It  is 
worth  noting  that  in  the  corresponding  experiments  with  small  LWC,  Figure  6.17 
showed  less  resemblance  to  Figure  6.15.  The  only  difference  between  Figures  6.22 
and  6.17  is  that  there  is  more  cloud  LWC  to  support  the  Bergeron-Findeison’s 
processes  in  the  former  than  in  the  latter  case.  As  a  result,  more  cloud  LWC  in 
Figure  6.22  is  converted  to  cloud  IWC,  which  precipitates  onto  the  ground  by  virtue 
of  gravitational  settling.  Consequently,  total  cloud  water  is  reduced  wherever  the 
cloud  LWC  is  large  enough  to  support  strong  Bergeron-Findeison’s  processes.  This 
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Figure  6.20.  Zonally  averaged  differences  in  cloud  LWC  with  a  small  autoconversion 
rate  (Co  =  10-,}sec-1). 
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Figure  6.22.  Zonally  averaged  differences  in  total  cloud  water  with  a  small  auto¬ 
conversion  rate  (Co  =  10"'4sec_1). 


may  explain  the  increase  of  total  cloud  water  in  Figure  6.17  due  to  the  small  cloud 
LWC  in  the  simulation. 

Comparisons  between  Figures  6.22  and  6.20  indicate  the  differences  in  the  total 
cloud  water  content  are  smaller  than  those  in  the  cloud  LWC,  since  the  existence 
of  cloud  ice  content  represents  part  of  the  depleted  cloud  LWC.  The  total  cloud 
water  content  decreases  in  the  region  where  ice  clouds  are  present.  The  reduction 
of  the  total  cloud  water  content  is  mainly  controlled  either  by  autoconversion  that 
converts  cloud  LWC  into  precipitation  or  by  gravitational  settling  that  removes  ice 
crystals  from  the  cloud  layer.  The  decrease  of  total  cloud  water  content  when  the 
ice  phase  processes  are  included  indicates  that  gravitational  settling  is  an  efficient 
mechanism  for  reducing  the  cloud  IWC. 
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Section  7 


SUMMARY  AND  CONCLUSION 

A  time-dependent,  three-dimensional,  large-scale  cloud-moisture  model  was  de¬ 
veloped  for  the  prediction  of  water  vapor,  cloud  LWC,  cloud  IWC,  precipitation  and 
temperature.  The  first  four  variables  are  governed  by  the  law  of  mass  conservation. 
The  latent  heat  exchanges  during  the  phase  transitions  of  the  water  substances 
are  ruled  by  the  first  law  of  thermodynamics  through  the  temperature  field.  The 
large-scale  cloud  cover  is  diagnostic  based  on  the  model-predicted  relative  humidity. 
Partial  cloudiness  is  allowed  to  form  when  large-scale  relative  humidity  is  less  than 
100%.  In  the  vertical,  a  16-layer  z-coordinate  is  adopted.  The  horizontal  resolution 
is  described  by  48  equally-spaced  longitudinal  points  and  38  Gaussian  latitudes. 

The  cloud  microphysical  processes  simulated  in  the  model  include  liquid  phase  as 
well  as  ice  phase,  which  is  physically  included  for  the  first  time  in  a  large-scale  cloud 
model  in  order  to  simulate  cirrus  clouds.  The  liquid  phase  processes  consist  of  evap¬ 
oration,  condensation,  autoconversion,  and  precipitation.  The  ice-phase  processes 
include  heterogeneous  nucleation  to  generate  ice  crystals,  depositional  growth  to 
simulate  Bergeron-Findeison’s  processes,  sublimation  to  deplete  ice  crystals,  and 
gravitational  settling  to  remove  ice  crystals  from  clouds.  The  Bergeron-Findeison’s 
processes  occur  primarily  at  temperatures  between  0°  C  and  —  40°  C,  while  homo¬ 
geneous  nucleation  dominates  below  —  40°  C.  Gravitational  settling  is  established 
by  mean  volume-weighted  downward  ice  flux,  which  is  obtained  by  integrating  the 
falling  velocity  and  mass  of  an  individual  ice  crystal  over  t.'ie  whole  spectrum  of  ice 
crystals. 
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The  radiative  transfer  program  is  based  on  a  broadband  method  and  involves 
the  transfer  of  IR  and  solar  radiation  in  clear  and  cloudy  regions.  In  a  cloudy 
atmosphere,  low  and  middle  clouds  are  treated  as  blackbodies  in  the  infrared 
radiative  transfer  calculations.  The  broadband  IR  emissivity,  reflectivity,  and 
transmissivity  for  cirrus  clouds  as  well  as  the  broadband  solar  absorption,  reflection, 
and  transmission  values  for  various  cloud  types  are  computed  based  on  the  param- 
eterizations  developed  by  Liou  and  Wittman(1979).  To  perform  radiative  transfer 
calculations  and  to  compare  with  observed  cloudiness,  the  model-predicted  cloud 
covers  are  strapped  into  low,  middle  and  high  clouds,  which  overlap  each  other  in 
a  statistically  random  manner. 

Large  amounts  of  satellite  data,  including  the  cloud  climatology  (3DNEPII) 
and  the  Earth  Radiation  Budget  (ERB)  data,  have  been  processed  into  formats 
appropriate  for  verification.  The  SMMR  data  consisting  of  vertically  integrated 
liquid  water  content  are  used  to  verify  model-predicted  cloud  LWC.  In  addition,  the 
cloud  cover  climatology  compiled  by  London  (1957)  based  on  surface  observations  is 
adopted  as  a  complementary  verification  dataset  to  assess  the  model-predicted  cloud 
cover  in  the  arctic  and  antarctic  regions,  where  the  satellite  observational  cloud 
covers  are  very  inaccurate  because  3DNEPH  has  not  been  able  to  distinguished 
cold  surfaces  from  clouds  aloft. 

A  96-hour  model  simulation,  with  the  initial  conditions  taken  at  12Z  on  July  1, 
1979,  has  been  carried  out  to  assess  the  performance  of  the  large-scale  cloud  model. 
The  wind  field  is  prescribed  using  the  results  from  the  AFGL  GCM.  Thus,  there  is 
no  cloud-radiation  feedback  to  the  dynamic  structure.  In  the  simulation,  the  cloud 
model  forecasts  cloud  cover,  cloud  LWC  and  IWC,  which  are  interactively  input 
into  the  radiation  model.  The  radiation  model  calculates  radiative  heating/cooling 
rates  based  on  the  prediction  cloud  properties.  The  radiative  heating/cooling  rates, 
in  turn,  feed  back  to  the  cloud  model  through  the  thermodynamic  equation,  and 
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modify  the  cloud  fields.  Thus,  cloud-radiation  interaction  through  the  medium  of 
the  cloud  cover,  cloud  LWC  and  IWC  is  accomplished  with  the  aid  of  the  large-scale 
cloud  model.  In  addition,  the  radiation  model  computes  OLR  and  net  solar  flux  at 
TOA  to  be  compared  with  the  ERB  products. 

The  comparisons  between  the  model  simulation  and  the  satellite  observations 
indicate  the  large-scale  cloud  model  is  capable  of  realistically  simulating  zonal 
means  and  geographical  distributions  of  cloud  fields,  including  cloud  cover,  OLR 
and  net  solar  flux  at  TOA,  and  cloud  LWC.  Differences  in  the  zonal  mean  total 
cloud  cover  between  the  simulation  and  the  3DNEPII  data  are  generally  within 
10%  of  total  cloudiness.  Generally,  the  simulated  OLR  and  net  solar  flux  at  TOA 
are  within  10  and  20  W/m2  of  the  ERB  results,  respectively.  The  zonal  mean 
values  of  simulated  cloud  LWP  compare  quite  well  with  the  SMMR  observations, 
with  differences  being  less  than  5  mg/cm2,  which  is  the  accuracy  of  the  SMMR  data. 
However,  it  is  clearly  indicated  that  the  model  underestimates  cloud  cover  and  cloud 
LWC  and  overestimates  the  OLR  in  the  tropics,  where  clouds  are  predominantly 
convective  clouds,  which  are  not  simulated  in  this  model. 

In  the  geographical  distributions  of  cloud  fields,  the  large-scale  cloud  model  has 
successfully  simulated  the  large  cloud  bands  in  the  ITCZ,  monsoon  areas,  SPCZ, 
and  the  mid-latitude  storm  tracks,  although  the  ITCZ  cloud  band  is  less  pronounced 
in  the  simulation.  At  the  subtropical  highs,  areas  of  small  cloudiness  over  southern 
Africa  and  northern  Africa,  and  North  and  South  America  are  clearly  identified  in 
both  of  the  observations  and  the  model  simulation.  The  simulated  precipitation  on 
the  east  coast  of  America  is  in  contrast  to  that  on  the  west  coast;  i.e.,  the  model 
produces  large  amounts  of  precipitation  on  the  east  coast  but  little  precipitation 
on  the  west  coast. 

The  model-simulated  cloud  fields  are  strongly  tied  to  large-scale  circulations. 

Strong  low-level  wind  convergences  and  deep  moist  tropospheres  in  the  ITCZ  and 
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monsoon  areas  assist  in  the  formation  of  large  cloud  amounts  in  these  regions. 
Great  baroclinic  instability  in  the  mid-latitude  storm  tracks  aids  in  the  formation 
of  large  amounts  of  cloud  in  these  areas.  The  clear-sky  conditions  at  the  subtropical 
highs  are  associated  with  the  large-scale  subsidence  in  the  branches  of  downward 
motion  in  the  Hardley  cell.  The  downward  motion  acts  to  prevent  water  vapor 
from  condensing  and  greatly  reduces  the  possibility  of  the  cloud  formation.  In 
brief,  upward  motion  supports  cloud  formation,  while  downward  motion  favors 
cloud  depletion. 

Observing  the  model-predicted  cloud  fields,  such  as  cloud  cover,  OLR  and  net 
solar  flux  at  TOA,  cloud  LWC  and  precipitation,  we  find  these  cloud  proper¬ 
ties  are  closely  related  one  another.  Large  amounts  of  cloud  cover  coincide  with 
abundant  cloud  LWC,  which  generates  large  precipitation  rates  by  virtue  of  the 
autoconversion  process  simulated  in  the  model.  In  addition,  the  model  simulation 
demonstrates  the  consistency  between  the  simulated  OLR  and  net  solar  flux  at 
TOA.  For  instance,  both  of  them  are  small  over  the  ITCZ  and  monsoon  areas, 
and  large  in  the  regions  of  the  subtropical  highs.  In  summary,  the  atmospheric 
moisture  and  cloud  fields  are  linked  together  by  the  cloud  microphysical  processes 
simulated  in  the  cloud  model.  More  importantly,  the  moisture,  cloud  fields  and 
earth’s  radiation  budget  are  physically  related  to  large-scale  thermal  and  dynamical 
structures,  such  as  temperature  and  vertical  velocity. 

The  magnitude  of  the  cloud  IWC  in  the  simulation  is  about  0.01  g/m3,  which 
is  consistent  with  the  aircraft  measurements  in  cirrus  clouds  as  well  as  with  other 
model  results.  The  maxima  of  cloud  IWC  occurs  along  the  —  15°  C  temperature 
line,  where  the  Bergeron-Findeison’s  processes  dominate.  This  result  agrees  well 
with  the  cloud-microphysics  laboratory  experiment  which  indicated  the  maximum 
ice  crystal  depositional  rate  occurred  at  —  15°  C  due  to  a  large  vapor  pressure 
difference  and  a  large  capacitance  value  at  this  temperature.  Because  of  the  strong 
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temperature  dependence  of  the  ice  crystal  growth  rate,  the  simulated  cloud  IWC 
appears  to  depend  more  on  temperature  than  on  vertical  velocity,  which  is  essential 
to  the  condensational  process. 

The  latitudinal  distribution  of  cloud  IW C  appears  to  be  realistic.  In  the  antarctic 
area,  all  of  the  low,  middle  and  high  clouds  are  ice  clouds.  However,  in  the  Arctic, 
only  middle  and  high  clouds  are  ice  clouds.  In  the  tropics,  only  high  clouds  contain 
ice  crystals.  The  simulated  cloud  LWCs  prevail  in  low  and  middle  clouds  at  the 
middle-low  latitudes  where  the  temperature  is  warm  and  abundant  moisture  is 
available  for  condensation. 

Sensitivity  studies  have  been  performed  to  examine  the  effects  of  radiative  heat¬ 
ing  fields  and  ice-phase  cloud  microphysics  on  large-scale  cloud  formation.  Numer¬ 
ical  experiments  have  been  carried  out  with  and  without  the  inclusion  of  radiation 
and  ice-phase  processes.  The  inclusion  of  radiative  processes  in  cloud  formation 
overwhelmingly  decreases  the  cloud  temperature,  with  the  strongest  cooling  oc¬ 
curring  at  the  middle  cloud  top.  However,  there  may  be  weak  heating  located 
in  a  low  cloud  if  cloud  cover  is  large  enough  to  produce  a  strong  downward  flux. 
The  decrease  of  temperature  by  radiative  cooling  increases  the  atmospheric  relative 
humidity  and  condensation  in  clouds.  As  a  result,  cloud  cover  and  total  cloud  water, 
including  cloud  LWC  and  IWC,  increase  due  to  the  inclusion  of  radiative  processes. 

The  inclusion  of  ice-phase  processes  has  a  greater  effect  on  the  cloud  LWC  and 
IWC  than  on  the  temperature  field.  The  differences  in  cloud  cover  due  to  the 
inclusion  of  ice  phase  are  insignificant.  Due  to  the  release  of  latent  heat  in  the 
Bergeron-Findeison’s  processes,  temperature  increases  occur  in  the  regions  where 
the  cloud  IWCs  prevail. 

The  changes  of  the  total  cloud  water  due  to  the  inclusion  of  the  ice  phase  is 
complicated  by  the  dependence  of  depositional  growth  on  the  amount  of  cloud 
LWC.  An  ice  crystal’s  growth  is  stopped  if  there  is  no  cloud  water  to  supply  water 
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vapor.  If  the  cloud  LWC  is  large  enough  to  support  effective  Bergeron-Fir.deison’s 
processes,  the  total  cloud  water  content  decreases  in  the  regions  where  ice  clouds 
dominate.  The  reduction  of  the  total  cloud  water  content  is  mainly  controlled  either 
by  autoconversion  which  converts  cloud  LWC  into  precipitation,  or  by  gravitational 
settling,  which  removes  ice  crystals  from  the  cloud  layer.  The  decrease  of  total  cloud 
water  content  when  the  ice-phase  processes  are  included  indicates  that  gravitational 
settling  is  an  efficient  mechanism  for  reducing  the  cloud  IWC. 
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APPENDIX  A 


DERIVATION  OF  PRECIPITATION  EQUATION  IN 

(2.5) 


The  mass  conservation  equation  for  rain  water,  qr ,  may  be  written  as  follow: 


c)  1  c/  d  d 


=  -prjP  -  (1  -  vi)?Er  .  (A.l) 

The  symbols  in  this  equation  are  defined  in  Chapter  2.  The  raindrop  terminal 
velocity  is  denoted  as  wr,  which  is  govern  by  the  gravitational  force  and  drag  force 
exerted  on  the  raindrop.  This  equation  may  be  further  simplified  based  on  the 
following  scale  arguments  consideration. 

The  magnitude  of  the  raindrop  terminal  velocity,  tty,  may  be  estimated  by 
assuming  a  balance  between  the  gravitational  force  and  drag  force  exerted  on  the 
raindrop.  Rogers  (1979)  derived  the  following  terminal  velocity  for  raindrops. 

=  2.2  x  103(po//3)1/2(r)1/2  ,  (A.2) 

where  p  is  the  air  density  and  po  is  the  reference  air  density  of  1.2  x  10-3  g/cm3, 
and  r  the  radius  of  a  raindrop.  At  700  mb,  p  is  roughly  estimated  as  0.9  x  10“3 
g/cm3.  Based  on  Equation  (A.2),  the  terminal  velocity  for  a  raindrop  with  0.6  mm 
radius  is  about  6  m/s,  i.e.,  on  the  order  of  10°  m/s. 

Scale  analysis  is  a  suitable  method  for  evaluating  the  magnitudes  of  various  terms 

in  the  governing  equations  for  a  particular  space-time  domain.  In  the  scale  analysis, 
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a  large-scale  motion  is  defined  by  the  following  characteristic  length,  depth,  and 
time  scales. 


L  ~  106  m  , 
V  ~  101  ms-1  , 
D  ~  104  m  , 
W  ~  10"2  ms"1  , 
T  =  L/V  ~  105  s"1  , 
Wr  ~  10°  ms-1  , 


where  L  and  D  denote  the  horizontal  and  vertical  geometry  scale,  respectively.  V, 
W,  and  T  are  the  typical  scale  values  of  horizontal  velocity,  vertical  velocity,  and  the 
life  span  of  large-scale  motions,  and  Wr  is  the  order  of  raindrop  terminal  velocity. 

We  can  now  estimate  the  magnitudes  of  the  following  terms  in  Equation  (A.l) 
as  follows: 

d  d 

~  Wr)}  ~  ^[pqr(-Wr)] 

^|J:«  10-v* 

~[(  Kl-Wr/D-lO-V1  . 

The  above  scale  analysis  shows  that  the  time  tendency  and  horizontal  advection 
terms  in  Equation  (A.l)  are  about  an  order  of  magnitude  smaller  than  the  down¬ 
ward  flux  of  precipitated  rain  water,  i.e.,  —pqrwr.  To  a  good  approximation,  the 
horizontal  advection  and  time  tendency  of  rain  water  may  be  neglected  due  to  the 
large  falling  velocity  of  rain  water.  Therefore,  the  rain  water  budget  equation  may 
be  rewritten  as  follows: 

Q 

~x~{P)  =  pr)P  -  (1  -  r))pET  ;  where  P  =  pqrwr. 
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APPENDIX  B 


DERIVATION  OF  LARGE-SCALE  CONDENSATION 


The  large-scale  condensation,  Q ,  is  related  to  the  time  rate  of  change  of  the 
saturation  vapor  pressure  as  follows: 

dqa 


Q  = 


dt 


where  the  saturation  specific  humidity,  q3,  is  given  by 

UT)  =  e-^P-  , 


(B.l) 


here  £  =  0.622  and  es(T )  is  a  saturation  vapor  pressure,  which  is  related  to  a 
large-scale  temperature  by  the  Clausius- Clapeyron  equation 

(B.2) 


dp 

—  =  LdT/RvT 2 

e. 


Differentiation  of  Equation  (B.l)  and  subsitution  from  Equation  (B.2)  leads  to 

1  dq,  L  dT  1  dP 


q3  dt  RVT2  dt  P  dt 
Furthermore,  the  thermodynamic  equation  gives 

_  dg, 
dt 


CpdT  RaT  dP 


(B.3) 


(B.4) 


L  dt  L  P  dt 
Eliminating  dT/dt  between  Equations  (B.3)  and  (B.4),  we  obtain 

dq1=-  *'LRa-CpRvT  u 
dt  Qs  [CPRVT2  +  qsL2  P  ’ 

where  u>  is  the  vertical  velocity  in  the  pressure  coordinate,  i.e.,  oj  =  dE-.  In 
large  scale  motions,  it  is  assumed  u>  ~  —pgw,  based  on  the  hydrostatic  equation. 
Therefore,  the  large-scale  condensation  in  our  model  is  given  by 
^  d  qs  n  q,{f)  r  LRa  -  CPRVT  ,  _ 

?=-¥=r'-RrW+5.i’1"  ’ 

where  Tj  is  the  dry  adiabatic  lapse  rate,  i.e.,  g/Cp. 
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APPENDIX  C 


PARAMETERIZATION  OF  RAINDROPS 
EVAPORATION 

For  a  unit  volume  of  air,  the  rate  of  evaporation  of  raindrops  can  be  obtained 
by  integrating  the  evaporation  of  an  individual  raindrop  over  the  whole  spectrum 
of  raindrops.  This  is  expressed  by 

(G1) 

where  /,  mr,  n  and  D  are  the  ventilation  factor,  the  mass  of  the  individual  raindrop, 
the  number  of  raindrops  per  unit  volume  of  air  per  unit  diameter  range  6D,  and 
the  raindrop  diameter,  respectively. 

The  rate  of  change  in  the  mass  of  an  individual  spherical  raindrop  may  be  derived 
from  Fick’s  law  and  formulated  as  follows: 

~  =  ItpDDM,  -  *0  ,  (C.2) 

where  Dd  is  the  mass  diffusivity  of  water  vapor  in  air. 

It  is  assumed  that  the  distribution  of  raindrops,  n,  follows  the  Marshall-Palmer 
(1948)  distribution  throughout  the  whole  process  of  evaporation.  This  distribution 
is  given  by 

n  =  nQ  e~XD  ,  (C.3) 

where  no  and  A  are  constants  which  are  determined  empirically.  Subsituting  Equa¬ 
tions  (C.2)  and  (C.3)  into  Equation  (C.l),  we  obtain 

Er  =  2nDd(qa  -  qQ)n0\~2  f  . 
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(C.4) 


To  relate  Er  with  the  total  precipitation  flux,  P ,  we  adopt  the  following  relation: 


P=  r wT^DzpwnQe~wdD  ,  (C.5) 

Jo  o 

where  wT  is  the  raindrop  terminal  velocity,  which  can  be  approximated  by  the 
following  equation  (Rogers,  1977), 

wr  =  2.2  x  103(/>o/p)1/2(£/2)1/2  ,  (C.6) 

where  p  is  the  air  density  and  p0  is  the  reference  air  density  of  1.2  x  10“3  g/cm3. 
Using  the  Marshall- Palmer  size  distribution,  we  obtain  the  following: 

A  =  4.529P-0'21  .  (C.7) 

Subsituting  Equation  (C.7)  into  Equation  (C.4),  leads  to 

Er  =  5.54  x  I0~3qs(l  —  ro)POA2  f  , 

where  we  have  used  a  value  of  0.226  cm2/sec,  for  the  water  vapor  diffusion  in  air, 
A/,  at  0°  C.  Let  Kg  =  5.54  x  10 ~3qsf]  then  we  have 

£r  =  /rE(l-r„)P°-4!  . 
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